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Preface. 

The present publication is an essay that was sent in (December, 1887) to com- 
pete for the Smith's Prizes at Cambridge. 

To the onlooker it is always a mournful thing to see what he considers splendid 
abilities or opportunities wasted for lack of knowledge of some paltry common- 
place truth. Such is in the main my feeling when considering the neglect of the 
study of Quaternions by that wonderful corporation the University of Cambridge. 
To the alumnus she is apt to appear as the leader in all branches of Mathematics. 
To the outsider she appears rather as the leader in Applied Mathematics and as a 
ready welcomer of other branches. 

If Quaternions were simply a branch of Pure Mathematics we could under- 
stand why the study of them was almost confined to the University which gave 
birth to them, but as the truth is quite otherwise it is hard to shew good reason 
why they have not struck root also in Cambridge. The prophet on whom Hamil- 
ton's mantle has fallen is more than a mathematician and more than a natural 
philosopher — he is both, and it is to be noted also that he is a Cambridge man. 
He has preached in season and out of season (if that were possible) that Quater- 
nions are especially useful in Physical applications. Why then has his Alma Mater 
turned a deaf ear? I cannot believe that she is in her dotage and has lost her hear- 
ing. The problem is beyond me and I give it up. 

But I wish to add my little efforts to Prof. Tait's powerful advocacy to bring 
about another state of affairs. Cambridge is the prepared ground on which if 
anywhere the study of the Physical applications of Quaternions ought to flourish. 

When I sent in the essay I had a faint misgiving that perchance there was not 
a single man in Cambridge who could understand it without much labour — and 
yet it is a straightforward application of Hamilton's principles. I cannot say what 
transformation scene has taken place in the five years that have elapsed, but an 
encouraging fact is that one professor at any rate has been imported from Dublin. 

There is no lack in Cambridge of the cultivation of Quaternions as an algebra, 
but this cultivation is not Hamiltonian, though an evidence of the great fecundity 
of Hamilton's work. Hamilton looked upon Quaternions as a geometrical method, 
and it is in this respect that he has as yet failed to find worthy followers resident 
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in Cambridge. [The chapter contributed by Prof. Cayley to Prof. Tait's 3rd ed. of 
'Quaternions' deals with quite a different subject from the rest of the treatise, a 
subject that deserves a distinctive name, say, Cayleyan Quaternions.] 

I have delayed for a considerable time the present publication in order at the 
last if possible to make it more effective. I have waited till I could by a more 
striking example than any in the essay shew the immense utility of Quaternions 
in the regions in which I believe them to be especially powerful. This I believe 
has been done in the 'Phil. Trans.' 1892, p. 685. Certainly on two occasions 
copious extracts have been published, viz. in the P. R. S. E., 1890-1, p. 98, and 
in the 'Phil. Mag.' June 1892, p. 477, but the reasons are simple. The first 
was published after the subject of the 'Phil. Trans.' paper had been considered 
sufficiently to afford clear daylight ahead in that direction, and the second after 
that paper had actually been despatched for publication. 

At the time of writing the essay I possessed little more than faith in the po- 
tentiality of Quaternions, and I felt that something more than faith was needed to 
convince scientists. It was thought that rather than publish in driblets it were bet- 
ter to wait for a more copious shower on the principle that a well-directed heavy 
blow is more effective than a long-continued series of little pushes. 

Perhaps more harm has been done than by any other cause to the study of 
Quaternions in their Physical applications by a silly superstition with which the 
nurses of Cambridge are wont to frighten their too timorous charges. This is the 
belief that the subject of Quaternions is difficult. It is difficult in one sense and 
in no other, in that sense in which the subject of analytical conies is difficult to 
the schoolboy, in the sense in which every subject is difficult whose fundamental 
ideas and methods are different from any the student has hitherto been introduced 
to. The only way to convince the nurses that Quaternions form a healthy diet 
for the young mathematician is to prove to them that they will "pay" in the first 
part of the Tripos. Of course this is an impossible task while the only questions 
set in the Tripos on the subject are in the second part and average one in two 
years. [This solitary biennial question is rarely if ever anything but an exercise 
in algebra. The very form in which candidates are invited, or at any rate were 
in my day, to study Quaternions is an insult to the memory of Hamilton. The 
monstrosity "Quaternions and other non-commutative algebras" can only be par- 
allelled by "Cartesian Geometry and other commutative algebras." When I was 
in Cambridge it was currently reported that if an answer to a Mathematical Tripos 
question were couched in Hebrew the candidate would or would not get credit for 
the answer according as one or more of the examiners did or did not understand 
Hebrew, and that in this respect Hebrew or Quaternions were strictly analogous.] 
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Is it hopeless to appeal to the charges? I will try. Let me suppose that some 
budding Cambridge Mathematician has followed me so far. I now address myself 
to him. Have you ever felt a joy in Mathematics? Probably you have, but it was 
before your schoolmasters had found you out and resolved to fashion you into an 
examinee. Even now you occasionally have feelings like the dimly remembered 
ones. Now and then you forget that you are nerving yourself for that Juggernaut 
the Tripos. Let me implore you as though your soul's salvation depended on it to 
let these trances run their utmost course in spite of solemn warnings from your 
nurse. You will in time be rewarded by a soul-thrilling dream whose subject is the 
Universe and whose organ to look upon the Universe withal is the sense called 
Quaternions. Steep yourself in the delirious pleasures. When you wake you will 
have forgotten the Tripos and in the fulness of time will develop into a financial 
wreck, but in possession of the memory of that heaven-sent dream you will be a 
far happier and richer man than the millionest millionaire. 

To pass to earth — from the few papers I have published it will be evident that 
the subject treated of here is one I have very much at heart, and I think that the 
publication of the essay is likely to conduce to an acceptance of the view that it 
is now the duty of mathematical physicists to study Quaternions seriously. I have 
been told by more than one of the few who have read some of my papers that they 
prove rather stiff reading. The reasons for this are not in the papers I believe but 
in matters which have already been indicated. Now the present essay reproduces 
the order in which the subject was developed in my own mind. The less complete 
treatment of a subject, especially if more diffuse, is often easier to follow than the 
finished product. It is therefore probable that the present essay is likely to prove 
more easy reading than my other papers. 

Moreover I wish it to be studied by a class of readers who are not in the 
habit of consulting the proceedings, &c, of learned societies. I want the slaves of 
examination to be arrested and to read, for it is apparently to the rising generation 
that we must look to wipe off the blot from the escutcheon of Cambridge. 

And now as to the essay itself. But one real alteration has been made. A pas- 
sage has been suppressed in which were made uncomplimentary remarks con- 
cerning a certain author for what the writer regards as his abuse of Quaternion 
methods. The author in question would no doubt have been perfectly well able 
to take care of himself, so that perhaps there was no very good reason for sup- 
pressing the passage as it still represents my convictions, but I did not want a side 
issue to be raised that would serve to distract attention from the main one. To 
bring the notation into harmony with my later papers dv and V which occur in 
the manuscript have been changed throughout to <i£ and A respectively. To fa- 
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cilitate printing the solidus has been freely introduced and the vinculum abjured. 
Mere slips of the pen have been corrected. A formal prefatory note required 
by the conditions of competition has been omitted. The Table of Contents was 
not prefixed to the original essay. It consists of little more than a collection of the 
headings scattered through the essay. Several notes have been added, all indicated 
by square brackets and the date (1892 or 1893). Otherwise the essay remains ab- 
solutely unaltered. The name originally given to the essay is at the head of p. 1 
below. The name on the title-page is adopted to prevent confusion of the essay 
with the 'Phil. Mag.', paper referred to above. What in the peculiar calligraphy of 
the manuscript was meant for the familiar fffO dg has been consistently rendered 
by the printer as fffO ds. As the mental operation of substituting the former for 
the latter is not laborious I have not thought it necessary to make the requisite 
extensive alterations in the proofs. 

I wish here to express my great indebtedness to Prof. Tait, not only for having 
through his published works given me such knowledge of Quaternions as I pos- 
sess but for giving me private encouragement at a time I sorely needed it. There 
was a time when I felt tempted to throw my convictions to the winds and follow 
the line of least resistance. To break down the solid and well-nigh universal scep- 
ticism as to the utility of Quaternions in Physics seemed too much like casting 
one's pearls — at least like crying in the wilderness. 

But though I recognise that I am fighting under Prof. Tait's banner, yet, as 
every subaltern could have conducted a campaign better than his general, so in 
some details I feel compelled to differ from Professor Tait. Some two or three 
years ago he was good enough to read the present essay. He somewhat severely 
criticised certain points but did not convince me on all. 

Among other things he pointed out that I sprung on the unsuspicious reader 
without due warning and explanation what may be considered as a peculiarity in 
symbolisation. I take this opportunity therefore of remedying the omission. In 
Quaternions on account of the non-commutative nature of multiplication we have 
not the same unlimited choice of order of the terms in a product as we have in 
ordinary algebra, and the same is true of certain quaternion operators. It is thus in- 
convenient in many cases to use the familiar method of indicating the connection 
between an operator and its operand by placing the former immediately before 
the latter. Another method is adopted. With this other method the operator may 
be separated from the operand, but it seems that there has been a tacit convention 
among users of this method that the separated operator is still to be restricted to 
precedence of the operand. There is of course nothing in the nature of things why 
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this should be so, though its violation may seem a trifle strange at first, just as the 
tyro in Latin is puzzled by the unexpected corners of a sentence in which adjec- 
tives (operators) and their nouns (operands) turn up. Indeed a Roman may be said 
to have anticipated in every detail the method of indicating the connection now 
under discussion, for he did so by the similarity of the suffixes of his operators 
and operands. In this essay his example is followed and therefore no restrictions 
except such as result from the genius of the language (the laws of Quaternions) 
are placed on the relative positions in a product of operators and operands. With 
this warning the reader ought to find no difficulty. 

One of Prof. Tait's criticisms already alluded to appears in the third edition 
of his 'Quaternions.' The process held up in § 500 of this edition as an exam- 
ple of "how not to do it" is contained in § 6 below and was first given in the 
'Mess, of Math.,' 1884. He implies that the process is a "most intensely artificial 
application of" Quaternions. If this were true I should consider it a perfectly le- 
gitimate criticism, but I hold that it is the exact reverse of the truth. In the course 
of Physical investigations certain volume integrals are found to be capable of, or 
by general considerations are obviously capable of transformation into surface 
integrals. We are led to seek for the correct expression in the latter form. Start- 
ing from this we can by a long, and in my opinion, tedious process arrive at the 
most general type of volume integral which is capable of transformation into a 
surface integral. [I may remark in passing that Prof. Tait did not however arrive 
at quite the most general type.] Does it follow that this is the most natural course 
of procedure? Certainly not, as I think. It would be the most natural course for 
the empiricist, but not for the scientist. When he has been introduced to one 
or two volume integrals capable of the transformation the natural course of the 
mathematician is to ask himself what is the most general volume integral of the 
kind. By quite elementary considerations he sees that while only such volume 
integrals as satisfy certain conditions are transformable into surface integrals, yet 
any surface integral which is continuous and applies to the complete boundary of 
any finite volume can be expressed as a volume integral throughout that volume. 
He is thus led to start from the surface integral and deduces by the briefest of 
processes the most general volume integral of the type required. Needless to say, 
when giving his demonstration he does not bare his soul in this way. He thinks 
rightly that any mathematician can at once divine the exact road he has followed. 
Where is the artificiality? 

Let me in conclusion say that even now I scarcely dare state what I believe to 
be the proper place of Quaternions in a Physical education, for fear my statements 
be regarded as the uninspired babblings of a misdirected enthusiast, but I cannot 
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refrain from saying that I look forward to the time when Quaternions will appear 
in every Physical text-book that assumes the knowledge of (say) elementary plane 
trigonometry. 

I am much indebted to Mr G. H. A. Wilson of Clare College, Cambridge, for 
helping me in the revision of the proofs, and take this opportunity of thanking 
him for the time and trouble he has devoted to the work. 

ALEX. MCAULAY. 

University of Tasmania, 

HOBART. 

March 26, 1893. 
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QUATERNIONS AS A PRACTICAL INSTRUMENT 
OF PHYSICAL RESEARCH. 

Section I. 
Introduction. 

It is a curious phenomenon in the History of Mathematics that the greatest 
work of the greatest Mathematician of the century which prides itself upon be- 
ing the most enlightened the world has yet seen, has suffered the most chilling 
neglect. 

The cause of this is not at first sight obvious. We have here little to do with 
the benefit provided by Quaternions to Pure Mathematics. The reason for the 
neglect here may be that Hamilton himself has developed the Science to such 
an extent as to make successors an impossibility. One cannot however resist a 
strong suspicion that were the subject even studied we should hear more from 
Pure Mathematicians, of Hamilton's valuable results. This reason at any rate 
cannot be assigned for the neglect of the Physical side of Quaternions. Hamilton 
has done but little in this field, and yet when we ask what Mathematical Physicists 
have been tempted by the bait to win easy laurels (to put the incentive on no higher 
grounds), the answer must be scarcely one. Prof. Tait is the grand exception to 
this. But well-known Physicist though he be, his fellow-workers for the most part 
render themselves incapable of appreciating his valuable services by studying 
the subject if at all only as dilettanti. The number who read a small amount 
in Quaternions is by no means small, but those who get further than what is 
recommended by Maxwell as imperatively necessary are but a small percentage 
of the whole. 

I cannot help thinking that this state of affairs is owing chiefly to a preju- 
dice. This prejudice is well seen in Maxwell's well-known statement — "I am 
convinced that the introduction of the ideas, as distinguished from the operations 
and methods of Quaternions, will be of great use to us in all parts of our sub- 
ject."* Now what I hold and what the main object of this essay is to prove is that 



* Elect, and Mag. Vol. I. § 10. 
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the "operations and methods" of Quaternions are as much better qualified to deal 
with Physics than the ordinary ones as are the "ideas". 

But, what has produced this notion, that the subject of Quaternions is only a 
pretty toy that has nothing to do with the serious work of practical Physics? It 
must be the fact that it has hitherto produced few results that appeal strongly to 
Physicists. This I acknowledge, but that the deduction is correct I strongly dis- 
believe. As well might an instrument of which nobody has attempted to master 
the principles be blamed for not being of much use. Workers naturally find them- 
selves while still inexperienced in the use of Quaternions incapable of clearly 
thinking through them and of making them do the work of Cartesian Geometry, 
and they conclude that Quaternions do not provide suitable treatment for what 
they have in hand. The fact is that the subject requires a slight development in 
order readily to apply to the practical consideration of most physical subjects. 
The first steps of this, which consist chiefly in the invention of new symbols of 
operation and a slight examination of their chief properties, I have endeavoured 
to give in the following pages. 

I may now state what I hold to be the mission of Quaternions to Physics. I 
believe that Physics would advance with both more rapid and surer strides were 
Quaternions introduced to serious study to the almost total exclusion of Carte- 
sian Geometry, except in an insignificant way as a particular case of the former. 
All the geometrical processes occurring in Physical theories and general Physical 
problems are much more graceful in their Quaternion than in their Cartesian garb. 
To illustrate what is here meant by "theory" and "general problem" let us take the 
case of Elasticity treated below. That by the methods advocated not only are the 
already well-known results of the general theory of Elasticity better proved, but 
more general results are obtained, will I think be acknowledged after a perusal 
of § 12 to §21 below. That Quaternions are superior to Cartesian Geometry in 
considering the general problems of ( 1 ) an infinite isotropic solid, (2) the torsion 
and bending of prisms and cylinders, (3) the general theory of wires, I have en- 
deavoured to shew in § 22-§ 33. But fox particular problems such as the torsion 
problem for a cylinder of given shape, we require of course the various theories 
specially constructed for the solution of particular problems such as Fourier's the- 
ories, complex variables, spherical harmonics, &c. It will thus be seen that I do 
not propose to banish these theories but merely Cartesian Geometry. 

So mistaken are the common notions concerning the pretensions of advo- 
cates of Quaternions that I was asked by one well-known Mathematician whether 
Quaternions furnished methods for the solution of differential equations, as he as- 
serted that this was all that remained for Mathematics in the domain of Physics! 
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Quaternions can no more solve differential equations than Cartesian Geometry, 
but the solution of such equations can be performed as readily, in fact generally 
more so, in the Quaternion shape as in the Cartesian. But that the sole work 
of Physical Mathematics to-day is the solution of differential equations I beg to 
question. There are many and important Physical questions and extensions of 
Physical theories that have little or nothing to do with such solutions. As witness 
I may call attention to the new Physical work which occurs below. 

If only on account of the extreme simplicity of Quaternion notation, large 
advances in the parts of Physics now indicated, are to be expected. Expressions 
which are far too cumbrous to be of much use in the Cartesian shape become so 
simple when translated into Quaternions, that they admit of easy interpretation 
and, what is perhaps of more importance, of easy manipulation. Compare for 
instance the particular case of equation (15m) § 16 below when F = with the 
same thing as considered in Thomson and Tait's Nat. Phil, App. C. The Quater- 
nion equation is 

PiSVupQwA = 0. 

The Cartesian exact equivalent consists of Thomson and Tait's equations (7), 
viz. 

d { _dw Ida ,\ dw da dw da 



2 — — + 1 + + 

dx y dA \dx / db dz dc dy 



+ 



d { dw da dw da dw (da 

— <2 + + — — + 1 

dy y dB dy da dz dc \ dx 



d ( dw da dw da dw [da 

+ —\2 + + — — + 1^ = 0, 

dz y dC dz da dy db \ dx 

and two similar equations. 

Many of the equations indeed in the part of the essay where this occurs, al- 
though quite simple enough to be thoroughly useful in their present form, lead 
to much more complicated equations than those just given when translated into 
Cartesian notation. 

It will thus be seen that there are two statements to make good: — (1) that 
Quaternions are in such a stage of development as already to justify the practically 
complete banishment of Cartesian Geometry from Physical questions of a general 
nature, and (2) that Quaternions will in Physics produce many new results that 
cannot be produced by the rival and older theory. 

To establish completely the first of these propositions it would be necessary 
to go over all the ground covered by Mathematical Physical Theories, by means 
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of our present subject, and compare the proofs with the ordinary ones. This of 
course is impossible in an essay. It would require a treatise of no small dimen- 
sions. But the principle can be followed to a small extent. I have therefore taken 
three typical theories and applied Quaternions to most of the general propositions 
in them. The three subjects are those of Elastic Solids, with the thermodynamic 
considerations necessary, Electricity and Magnetism, and Hydrodynamics. It is 
impossible without greatly exceeding due limits of space to consider in addition, 
Conduction of Heat, Acoustics, Physical Optics, and the Kinetic Theory of Gases. 
With the exception of the first of these subjects I do not profess even to have at- 
tempted hitherto the desired applications, but one would seem almost justified 
in arguing that, since Quaternions have been found so applicable to the subjects 
considered, they are very likely to prove useful to about the same extent in similar 
theories. Again, only in one of the subjects chosen, viz., Hydrodynamics, have 
I given the whole of the general theory which usually appears in text-books. For 
instance, in Electricity and Magnetism I have not considered Electric Conduc- 
tion in three dimensions which, as Maxwell remarks, lends itself very readily to 
Quaternion treatment, nor Magnetic Induction, nor the Electro-Magnetic Theory 
of Light. Again, I have left out the consideration of Poynting's theories of Elec- 
tricity which are very beautifully treated by Quaternions, and I felt much tempted 
to introduce some considerations in connection with the Molecular Current the- 
ory of Magnetism. With similar reluctance I have been compelled to omit many 
applications in the Theory of Elastic Solids, but the already too large size of the 
essay admitted of no other course. Notwithstanding these omissions, I think that 
what I have done in this part will go far to bear out the truth of the first proposition 
I have stated above. 

But it is the second that I would especially lay stress upon. In the first it is 
merely stated that Cartesian Geometry is an antiquated machine that ought to be 
thrown aside to make room for modern improvements. But the second asserts 
that the improved machinery will not only do the work of the old better, but 
will also do much work that the old is quite incapable of doing at all. Should 
this be satisfactorily established and should Physicists in that case still refuse 
to have anything to do with Quaternions, they would place themselves in the 
position of the traditional workmen who so strongly objected to the introduction 
of machinery to supplant manual labour. 

But in a few months and synchronously with the work I have already de- 
scribed, to arrive at a large number of new results is too much to expect even 
from such a subject as that now under discussion. There are however some few 
such results to shew. I have endeavoured to advance each of the theories chosen 
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in at least one direction. In the subject of Elastic Solids I have expressed the 
stress in terms of the strain in the most general case, i.e. where the strain is not 
small, where the ordinary assumption of no stress-couple is not made and where 
no assumption is made as to homogeneity, isotropy, &c. I have also obtained the 
equations of motion when there is given an external force and couple per unit 
volume of the unstrained solid. These two problems, as will be seen, are by no 
means identical. In Electrostatics I have considered the most general mechani- 
cal results flowing from Maxwell's theory, and their explanation by stress in the 
dielectric. These results are not known, as might be inferred from this mode of 
statement, for to solve the problem we require to know forty-two independent 
constants to express the properties of the dielectric at a given state of strain at 
each point. These are the six coefficients of specific inductive capacity and their 
thirty-six differential coefficients with regard to the six coordinates of pure strain. 
But, as far as I am aware, only such particular cases of this have already been 
considered as make the forty-two constants reduce at most to three. In Hydrody- 
namics I have endeavoured to deduce certain general phenomena which would be 
exhibited by vortex-atoms acting upon one another. This has been done by exam- 
ination of an equation which has not, I believe, been hitherto given. The result of 
this part of the essay is to lead to a presumption against Sir William Thomson's 
Vortex- Atom Theory and in favour of Hicks 's. 

As one of the objects of this introduction is to give a bird's-eye view of the 
merits of Quaternions as opposed to Cartesian Geometry, it will not be out of 
place to give side by side the Quaternion and the Cartesian forms of most of the 
new results I have been speaking about. It must be premised, as already hinted, 
that the usefulness of these results must be judged not by the Cartesian but by the 
Quaternion form. 



Elasticity. 

Let the point (x,y,z) of an elastic solid be displaced to (x',y',z'). The strain 
at any point that is caused may be supposed due to a pure strain followed by a 
rotation. In Section III. below, this pure strain is called if/. Let its coordinates be 
e, f, g, a/2, b/2, c/2; i.e. if the vector (t;,ri,D becomes (g,rf,g) by means of 
the pure strain, then 

&c, &c. 
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Thus when the strain is small e, f, g reduce to Thomson and Tait's 1 + e, 1 + /, 
1 + g and a, b, c are the same both in their case and the present one. Now let 
the coordinates of W, § 16 below, be E, F, G, A/2, 5/2, C/2. Equation (15), § 16 
below, viz. 

**¥oj = i]f 2 oj= X 'xu = ViSp\p' 2 Sa>V 2 , 
gives in our present notation 

E = e 2 + c 2 /4 + b 2 /4 = (dx'/dxf + (dy'/dxf + (dz/dxf, 

&c, &c. 
A = a(f + g) + Z>c/2 

= 2{{dx /dy){dx jdz) + (dy I dy)(dy' I dz) + (dz /dy)(dz /dz)}, 

&c, &c. 

which shew that the present £, F, G, A/2, 5/2, C/2 are the A, 6, C, a, b, c of 
Thomson and Tait's Nat. Phil., App. C. 
Let us put 

j{ xVz '\ = j 

\xyz 
j{ yZ \ = J n , &c, &c, 

J \) = Jn,J\ yZ ) = Jn, &c, &c, &c, &c. 

I have shewn in § 14 below that the stress-couple is quite independent of the 
strain. Thus we may consider the stress to consist of two parts — an ordinary 
stress PQRSTU as in Thomson and Tait's Nat. Phil, and a stress which causes a 
couple per unit volume L'M'N' . The former only of these will depend on strain. 
The result of the two will be to cause a force (as indeed can be seen from the 
expressions in § 13 below) per unit area on the x-interface P, U + N' /2, T - M' /2, 
and so for the other interfaces. If L, M, N be the external couple per unit volume 
of the unstrained solid we shall have 

L' = -L/J, M' = -M/ J, N' = -N/J, 



2 

This result is one of Tait's (Quaternions § 365 where he has </> <p = vj ). It is given here for 
completeness. 
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for the external couple and the stress-couple are always equal and opposite. Thus 
the force on the x-interface becomes 

P, U - N/2J, T + M/2J 

and similarly for the other interfaces. 

To express the part of the stress (P &c.) which depends on the strain in terms 
of that strain, consider w the potential energy per unit volume of the unstrained 
solid as a function of E &c. In the general thermodynamic case w may be defined 
by saying that 

w X (the element of volume) 

= (the intrinsic energy of the element) 

- (the entropy of the element X its absolute temperature X Joule's coefficient). 

Of course w may be, and indeed is in § 14, § 15 below, regarded as a function 
of e &c. 

The equation for stress is (15b) § 16 below, viz., 

Jipco = 2xyQwx'a> = 2p\S p' 2 (oSV HpQvvV \. 

The second of the expressions is in terms of the strain and the third in terms 
of the displacement and its derivatives. In our present notation this last is 

/ \ 2 i \2 / \ 2 

JP dx'\ dw dx'\ dw dx'\ dw 



2 \dx J dE \dy } dF \dz j dG 

„ dx' dx' dw „ dx' dx' dw „ dx' dx' dw 

+ 2 + 2 +2 , 

dy dz dA dz dx dB dx dy dC 

&c, &c. 

JS dy' dz' dw dy' dz' dw dy' dz' dw 

2 dx dx dE dy dy dF dz dz dG 

'dy' dz' dy' dz' \ dw (dy' dz' dy' dz' \ dw 
dy dz dz dy ) dA \ dz dx dx dz ) dB 

t dy' dz' dy' dz' \ dw 
\ dx dy dy dx ) dC ' 
&c, &c. 
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In § 14 I also obtain this part of the stress explicitly in terms of e, f, g, a, b, 
c, of w as a function of these quantities and of the axis and amount of rotation. 
But these results are so very complicated in their Cartesian shape that it is quite 
useless to give them. 

To put down the equations of motion let X x , Y x , Z x be the force due to stress 
on what before strain was unit area perpendicular to the axis of x. Similarly for 
X y , &c. Next suppose that X, Y, Z is the external force per unit volume of the 
unstrained solid and let D be the original density of the solid. Then the equation 
of motion (I5n) § 16a below, viz. 

Dp' = F + tA, 

gives in our present notation 

X + dX x /dx + dXy/dy + dX z /dz = x'D, &c, &c. 

It remains to express X x &c. in terms of the displacement and LMN. This is 
done in equation (15/) § 16 below, viz. 

to) = -2p\SVwawu) + 3VMVp' l p' 2 ScoViV 2 /2SViV 2 V 3 Sp' l p' 2 p' 3 : 
In our present notation this consists of the following nine equations: 



-A 



dw dx' dw dx' dw dx' \ J]jN - J^M 



\dE dx dC dy dB dz ) 2/ 

tdwdy' dw dy' dw dy'\ J U L - J U N 



Y x = 2\ — + - + — + 

\dE dx dC dy dB dz ) 2J 



Z, 



dw dz dw dz dw dz \ J\\M - J n L 

= 2 + + + , 

\dE dx dC dy dB dz j 27 



and six similar equations. 

We thus see that in the case where LMN are zero, our present X x , X y , X z are 
the PQR of Thomson and Tait's Nat. Phil. App. C (d), and therefore equations (7) 
of that article agree with our equations of motion when we put both the external 
force and the acceleration zero. 



The second term on the right contains in full the nine terms corresponding to (JnN - 
JnM)/2J. Quaternion notation is therefore here, as in nearly all cases which occur in Physics, 
considerably more compact even than the notations of determinants or Jacobians. 
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These are some of the new results in Elasticity, but, as I have hinted, there 
are others in § 14, § 15 which it would be waste of time to give in their Cartesian 
form. 

Electricity. 

In Section IV. below I have considered, as already stated, the most general 
mechanical results flowing from Maxwell's theory of Electrostatics. I have shewn 
that here, as in the particular cases considered by others, the forces, whether per 
unit volume or per unit surface, can be explained by a stress in the dielectric. It is 
easiest to describe these forces by means of the stress. 

Let the coordinates of the stress be PQRSTU. Then FiF 2 F 3 the mechanical 
force, due to the field per unit volume, exerted upon the dielectric where there is 
no discontinuity in the stress, is given by 

Fi = dP/dx + dU/dy + dT/dz, &c, &c. 

and (/, m, n) being the direction cosines of the normal to any surface, pointing 
away from the region considered 

F[ = -[IP + mU + nT] a -[] b , &c, &c, 

where a, b indicate the two sides of the surface and F{ , F 2 ', F 3 ' is the force due 
to the field per unit surface. 

It remains to find P &c. Let X, Y, Z be the electro-motive force, a, fi, y the 
displacement, w the potential energy per unit volume and K xx , K yy , K zz , K yz , K zx , 
K xy the coefficients of specific inductive capacity. Let 1 + e, 1 + /, 1 + g, a/2, 
b/2, c/2 denote the pure part of the strain of the medium. The ^'s will then be 
functions of e &c. and we must suppose these functions known, or at any rate we 
must assume the knowledge of both the values of the ^'s and their differential 
coefficients at the particular state of strain in which the medium is when under 
consideration. The relations between the above quantities are 

Ana = K XX X + K xy Y + K ZX Z, &c, &c. 
w = (Xa + Yfi + Zy)/2 

= (K XX X 2 + K yy Y 2 + K ZZ Z 2 + 2K yz YZ + 2K ZX ZX + 2K xy XY)/Sn. 

It is the second of these expressions for w which is assumed below, and the 
differentiations of course refer only to the ^'s. The equation expressing P &c. in 
terms of the field is (21) § 40 below, viz. 

<pco = -\V~DoM - ipQww, 
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which in our present notation gives the following six equations 

p = -\{-aX+pY + yZ)-dwlde, Sec, Sec, 
S = \(fiZ + yY) - dw/da, Sec, Sec 

I have shewn in § 41— § 45 below that these results agree with particular results 
obtained by others. 

Hydrodynamics. 

The new work in this subject is given in Section VI. — "The Vortex-Atom 
Theory." It is quite unnecessary to translate the various expressions there used 
into the Cartesian form. I give here only the principal equation in its two chief 
forms, equation (9) § 89 and equation (1 1) § 90, viz. 

P + v - a 2 12 + (Any 1 fff (ScttVu + udm/dt) rfs = H, 

P + v- a 2 12 + {Any 1 fff{dsS Vu(V(tt - ma) + ud(mds)/dt} = H. 

In Cartesian notation these are 



I 



dp/p + V + q 2 /2 

~ W'jOJPK*' - *X W? 7 - v£) + ••• + •• • ]/r 3 
+ (dc/dt)/r) dx' dy' dz = H. 

dp/p + V + q 2 /2 

~ ( 4n y l fff{( x ' - x)[2(wt] - vQ - cu] + ••• + •• • }/r 3 . dx' dy' dz 
~ (4ny l fff{d(cdx' dy' dz')/dt}/r = H. 



The fluid here considered is one whose motion is continuous from point to 
point and which extends to infinity. The volume integral extends throughout 
space. The notation is as usual. It is only necessary to say that H is a function of 
the time only, r is the distance between the points x' , y',z' and x, y, z; 

c = du/dx + dv/dy + dw/dz; 

d/dt is put for differentiation which follows a particle of the fluid, and d/dt for 
that which refers to a fixed point. 

The explanation of the unusual length of this essay, which I feel is called for, 
is contained in the foregoing description of its objects. If the objects be justifiable, 
so must also be the length which is a necessary outcome of those objects. 



Section II. 
Quaternion Theorems. 

Definitions. 

1. As there are two or three symbols and terms which will be in constant 
use in the following pages that are new or more general in their signification 
than is usual, it is necessary to be perhaps somewhat tediously minute in a few 
preliminary definitions and explanations. 

A function of a variable in the following essay is to be understood to mean 
anything which depends on the variable and which can be subjected to mathemat- 
ical operations; the variable itself being anything capable of being represented by 
a mathematical symbol. In Cartesian Geometry the variable is generally a sin- 
gle scalar. In Quaternions on the other hand a general quaternion variable is not 
infrequent, a variable which requires 4 scalars for its specification, and similarly 
for the function. In both, however, either the variable or the function may be a 
mere symbol of operation. In the following essay we shall frequently have to 
speak of variables and functions which are neither quaternions nor mere symbols 
of operation. For instance K in § 40 below requires 6 scalars to specify it, and it 
is a function of ip which requires 6 scalars and p which requires 3 scalars. When 
in future the expression "any function" is used it is always to be understood in the 
general sense just explained. 

We shall frequently have to deal with functions of many independent vectors, 
and especially with functions which are linear in each of the constituent vectors. 
These functions merely require to be noticed but not defined. 

Hamilton has defined the meaning of the symbolic vector V thus: — 

„ d d d 
dx dy dz 

where i, j, k are unit vectors in the directions of the mutually perpendicular axes x, 
y, z. I have found it necessary somewhat to expand the meaning of this symbol. 
When a numerical suffix 1, 2, ... is attached to a V in any expression it is to 
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indicate that the differentiations implied in the V are to refer to and only to other 
symbols in the same expression which have the same suffix. After the implied 
differentiations have been performed the suffixes are of course removed. Thus 
Q(a,/3,y,6) being a quaternion function of any four vectors a, /3, y, 5, linear in 
each 

„ „ IdAdu. \ IdAdu \ IdAdu 



dxdx J \dy dx j \ dz dx 

IdAdfi \ IdAdfi \ IdAdfi 

\dx dy ) \dy dy ) \dz dy 

IdAdfi \ IdAdu \ IdAdfi 

\dx dz j \dy dz J \dz dz 

and again 

Gi(^i,j"2,V 1 ,V 2 ) = G3Ui./*2,Vi + V 3 ,V 2 ). 

It is convenient to reserve the symbol A for a special meaning. It is to be 
regarded as a particular form of V, but its differentiations are to refer to all the 
variables in the term in which it appears. Thus Q being as before 

Q{AuH, A, VO = Q 2 (Aulu 2 , Vj + V 2 , Vj)*. 



If in a linear expression or function Vj and p! (p being as usual = ix + /y + fe) 
occur once each they can be interchanged. Similarly for V 2 andp 2 . So often does 
this occur that I have thought it advisable to use a separate symbol £,\ for each of 
the two Vi and p\, £ 2 for each of the two V 2 and p 2 and so for £3, &c. If only 
one such pair occur there is of course no need for the suffix attached to £. Thus £ 
may be looked upon as a symbolic vector or as a single term put down instead of 
three. For Q(a,/3) being linear in each of the vectors a,fi 

Q({,0 = G(Vi,pi) = Q(i,i) + Q(jJ) + Q(k,k). (1) 

There is one more extension of the meaning of V to be given, u, v, w being 
the rectangular coordinates of any vector cr, ^V is defined by the equation 

d d d 

^V = i— + j— + k—. 
du dv dw 



'These meanings for Vi, V2 . . . A I used in a paper on "Some General Theorems in Quaternion 
Integration," in the Mess, of Math. Vol. xiv. (1884), p. 26. The investigations there given are for 
the most part incorporated below. [Note added, 1892, see preface as to the alteration of V into 
A.] 
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To o-V of course are to be attached, when necessary, the suffixes above explained 
in connection with V. Moreover just as for Vi, p\ we may put £, £ so also for o-Vi, 
<Ti may we put the same. 

With these meanings one important result follows at once. The V\'s, V 2 's, 
&c, obey all the laws of ordinary vectors whether with regard to multiplication or 
addition, for the coordinates d/dx, d/dy, d/dz of any V obey with the coordinates 
of any vector or any other V all the laws of common algebra. 

Just as o-V may be defined as a symbolic vector whose coordinates are d/du, 
d/dv, d/dw so <p being a linear vector function of any vector whose coordinates 
are 

{a\h\C\ a 2 b2C2 a^b^c^) (i.e. <f>i = a\i + b\j + c\k, &c). 

^Q* is defined as a symbolic linear vector function whose coordinates are 

(d/dai, d/dbi, d/dc\, d/da 2 , djdb 2 , d/dc 2 , d/da^,, d/db?,, d/dc?,), 

and to ^Q is to be applied exactly the same system of suffixes as in the case of V. 
Thus q being any quaternion function of (p, and to any vector 

tf,(Jito .q\=- (idq/dai + jdq/dbi + kdq/dci)Sia> 

- (idq/da 2 + jdq/db 2 + kdq/dc 2 )S jco 

- {idqlda-i + jdq/db^ + kdq/dc^Sku). 

The same symbol ^(j is used without any inconvenience with a slightly differ- 
ent meaning. If the independent variable be a self-conjugate linear vector func- 
tion it has only six coordinates. If these are PQRSTU (i.e. (pi = Pi+Uj+Tk,&c.) 
^Q is defined as a self-conjugate linear vector function whose coordinates are 

(41 dP, d/dQ, d/dR, \d/dS, \d/dT, \d/dU). 

We shall frequently have to compare volume integrals with integrals taken 
over the bounding surface of the volume, and again surface integrals with inte- 
grals taken round the boundary of the surface. For this purpose we shall use the 
following notations for linear, surface and volume integrals respectively J Q dp, 
JjQ dL, JJjQ ds where Q is any function of the position of a point. Here dp is a 
vector element of the curve, dL a vector element of the surface, and ds an element 
of volume. When comparisons between line and surface integrals are made we 
take dL in such a direction that dp is in the direction of positive rotation round 
the element dL close to it. When comparisons between surface and volume inte- 
grals are made dL is always taken in the direction away from the volume which it 
bounds. 



*I have used an inverted D to indicate the analogy to Hamilton's inverted A. 
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Properties of '£. 

2. The property of £ on which nearly all its usefulness depends is that if <x 
be any vector 

o- = -£S£cr, 

which is given at once by equation (1) of last section. 

This gives a useful expression for the conjugate of a linear vector function of 
a vector. Let cp be the function and co, r any two vectors. Then cp' denoting as 
usual the conjugate of <p we have 

Scocpr = Srcp'co, 

whence putting on the left t = —£S£t we have 

St(-^ScoM) = Srcp'co, 

or since t is quite arbitrary 

cp' co = -£Scocp£- (2) 

From this we at once deduce expressions for the pure part cpco and the rota- 
tional part Veco of cpco by putting 

(0 + $)co = -ftS co£ - £S coM = ^co, | 

(cp - cb')co = -MS to£ + £S coM = VV£M-u = 2Veco. J 

And all the other well-known relations between cp and cp' are at once given e.g. 
S£M = S %cb'%, i-e. the "convergence" of cp = the "convergence" of cp'. 

3. Let Q(A,/u) be any function of two vectors which is linear in each. Then 
if cpco be any linear vector function of a vector co given by 

cpco = -ZfiScoa | 

we have Q(£,M) = XQ(-£S£a,/3) = XQ(a,/3), J ^ 

or more generally 

Q(£,<f>xO = ZQ(x'a,f3). (4a) 

To prove, it is only necessary to observe that 

<Px£=-Y.pSa X £ = -?<pS&<x, 
and that -£S i,x a = X a - 
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(5) 



As a particular case of eq. (4) let (p have the self-conjugate value 

(pa) = -^LifiScoa + aStofi). 
Then <2(<T, H) = \*{Q(a,p) + Q(fi, a)}, 

or if Q(A, n) is symmetrical in A and fi 

Q({,<PO = ZQ(a,P). (6) 

The application we shall frequently make of this is to the case when for a we 
put Vi and for f3, cr\, where cr is any vector function of the position of a point. In 
this case the first expression for (p is the strain function and the second expression 
the pure strain function resulting from a small displacement a at every point. As 
a simple particular case put Q(A,fi) = S A/u so that Q is symmetrical in A and //. 
Thus <p being either of these functions 

Another important equation is 

'Q({,<f>0 = Q(<f>'{,0. (6a) 

This is quite independent of the form of (p. To prove, observe that by equa- 
tion (2) 

and that -£S£(p%\ = <p£,\. Thus we get rid of £ and may now drop the suffix of £i 
and so get eq. (6a). [Notice that by means of (6a), (4a) may be deduced from (4); 
for by (6a) 

g(£, <PxO = Q(x'^ 4>0 = ZQ(x'a,/3) by (4).] 

3a. A more important result is the expression for (p~ x a> in terms of cp. We 
assume that 

S 4>A (p/u (pv = mS A/jv, 

where A, /u,v are any three vectors and m is a scalar independent of these vectors. 
Substituting £i, g 2 , £3 for A, H, v and multiplying by Sgifrb 

S{ 1 hkSMiM2M3 = 6m, (6b) 



*[Note added, 1892. For practice it is convenient to remember this in words: — A term in 
which £ and <p£ occur is unaltered in value by changing them into <p'£ and £ respectively.] 
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which gives m in terms of cp- That S£\ & &S£i £2 £3 = 6 is seen by getting rid of 
each pair of £'s in succession thus: — 

S£ (S£ V&6K26 = -SW 3 . 66 = S(& 3 - &&&)& = -2^3 = 6. 

Next observe that 

S cpco (p£\ 4>^2 = mS a> £1 £2- 
Multiplying by V(i& and again on the right getting rid of the £s we have 

V£i£ S cpco cp£\ cp£ 2 = -2mco, (6c) 

whence from equation (6b) 

coS 666 S <p£i 06 06 = -3 V66 S <f>a> 06 06, 
or changing co into _1 cj 

, 3766^0606 (t ,~ 

d> co = . (6a) 

S666S060606 
By equation (6a) of last section we can also put this in the form 

, 3y0'6 0'6^66 rf . 

d> co = , (be) 

s 666^6 ^6 ^6 

so that cp~ l co is obtained explicitly in terms of cp or cp'. 

Equation (6c) or (6c?) can be put in another useful form which is more anal- 
ogous to the ordinary cubic and can be easily deduced therefrom, or *less easily 
from (6d), viz. 

S666(0 3 wS666 ~ 34> 2 coS^ 2 <Ph 

+ 3chco S 6 06 06 ~ coS 06 06 06) = 0. 

As a useful particular case of equation (6c?) we may notice that by equation (4) 
§3 if 



<pco = -CTiStoVi, 

-1 

and (p" 1 co = -3Vo- l o- 2 ScoV l V 2 /SV l V2 , V 3 S 0-^20-3. (6h) 



Cp l CO = -3W1V2S CO O-1CT2/SV \V 2^3 So-iCT2CT 3 , (6g) 



*[Note added, 1892. The cubic may be obtained in a more useful form from the equation 
wShh&SMi Hi Hi = -3V(i( 2 S<fiQi<Ki Hi thus 

VM 2 Scl>toH Hi = <P«>S . HH1VM1 - HxS . ^HiVMi + HiS . <f>coHiV^ 2 
= Hos . h Hi vhb - 2H s . Ho Hi VMi- 
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4. Let (p, if/ be two linear vector functions of a vector. Then if 

s x SK = SxCW, 

where ^ is a quite arbitrary linear vector function 

(p = iff, 
for we may put;^ = tS£u> where t and oj are arbitrary vectors, so that 

STcpco = Srif/to, 
or (pco = if/to. 

Similarly* if (p and if/ are both self-conjugate and ^ is a quite arbitrary self- 
conjugate linear vector function the same relation holds as can be seen by putting 

X( = tS £w + cjS (t. 

Fundamental Property of (J. 

5. Just as the fundamental property of o-V is that, Q being any function of <x 

6Q = -Q 1 S6<r < p u 

so we have a similar property of Q. Q being any function of a linear vector 
function 

6Q = -QiS6MtQ 1 £. (7) 



Again fa S . </>w fa V£i& - fa s ■ 4>uV ■ fa Vfrfr 

= fa Sfa(faS& fa + { 2 Sd fa) 
= -m S& fa)S& fa + <K(i S(\ fa)S& fa 
= (p 2 ajS£<t>Z;-(p 2 £S£(pa> = <p 2 coS£<f>£ + <p 3 co. 
Hence <p a> — m"<p a> + m fa — ma> = 0, 

where 6ra = S £ £> (3 S fa fa fa 

2m' =-SV{faVfa fa 
m" = -S£fa] 



"[Note added, 1892. The following slightly more general statement is a practically much 
more convenient form of enunciation: if Sx£, <f>£ = SxC &£> where ^ is a perfectly arbitrary self- 
conjugate and (/> and if/ are not necessarily self-conjugate then (f> = if/]. 
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The property is proved in the same way as for V, viz. by expanding S 6<f>g$Qig 
in terms of the coordinates of ^,q. First let <p be not self-conjugate, and let its nine 
coordinates be 

(a\b\C\ a 2 b 2 c 2 GI3&3C3). 

Thus 

-<2i S <fy£*CIi£ = -Gi S S(pi^a\i -Q\S SQJiQiJ -Q\S ScpktQik, 
= 6a\ dQ/dai + 8b\ dQ/dbi + 8c\ dQ/dc\, 
+ 5a 2 dQ/da 2 + 5b 2 dQ/db 2 + 6c 2 dQ/dc 2 , 
+ 5a-i dQ/da 3 + 5b?, dQ/dbi + SC3 dQIdc? = 6Q. 

The proposition is exactly similarly proved when <p is self-conjugate. 

Theorems in Integration. 

6. Referring back to § 1 above for our notation for linear surface and volume 
integrals we will now prove that if Q be any linear function of a vector* 

fQdp = ffQ(VdlA), (8) 

ffQdZ = fffQAds. (9) 

To prove the first divide the surface up into a series of elementary parallelo- 
grams by two families of lines — one or more members of one family coinciding 
with the given boundary, — apply the line integral to the boundary of each paral- 
lelogram and sum for the whole. The result will be the linear integral given in 
equation (8). Let the sides of one such parallelogram taken in order in the positive 
direction be a, (3 +/?', -a - a', -/?; so that a' and/?' are infinitely small compared 
with a and f3, and we have the identical relation 

= a +J3 + /3' - a - a' -j3 =ft - a . 

The terms contributed to J Q dp by the sides a and -a -a' will be (neglecting 
terms of the third and higher orders of small quantities) 

Qa-Qa- Qa' + QiaSfiV^ = -Qa' + QkxS/Nl 



'These two propositions are generalisations of what Tait and Hicks have from time to time 
proved. They were first given in the present form by me in the article already referred to in § 1 
above. In that paper the necessary references are given. 



§7.] QUATERNION THEOREMS. 19 

Similarly the terms given by the other two sides will be 

Qfi'-Q 1 J3SaV l 

so that remembering that /?' - a' = and therefore Q/3' - Qa' = we have for 
the whole boundary of the parallelogram 

QiaSpV, - Gu85orVi = Q^VVafi . Vj) = QVdZA, 

where <TL is put for Vaj3. Adding for the whole surface we get equation (8). 

Equation (9) is proved in an exactly similar way by splitting the volume up 
into elementary parallelepipeda by three families of surfaces one or more mem- 
bers of one of the families coinciding with the given boundary. If a, fi, y be the 
vector edges of one such parallelepiped we get a term corresponding to Qfi - Qa' 
viz. 

<2(vector sum of surface of parallelepiped) = 0, 

and we get the sum of three terms corresponding to 

Q 1 aS/3V 1 -Q l fiSaV 1 
above, viz. 

-Q l (Vfiy)SaV 1 - Q 1 (Vya) S fiV \ - Q l (Va/3)SyV l = -Q^.Safiy, 
whence putting Saj3y = -<is we get equation (9). 

7. It will be observed that the above theorems have been proved only for 
cases where we can put dQ = -Q\S dpV \ i.e. when the space fluxes of Q are 
finite. If at any isolated point they are not finite this point must be shut off from 
the rest of the space by a small closed surface or curve as the case may be and 
this surface or curve must be reckoned as part of the boundary of the space. If 
at a surface (or curve) Q has a discontinuous value so that its derivatives are 
there infinite whereas on each side they are finite, this surface (or curve) must be 
considered as part of the boundary and each element of it will occur twice, i.e. 
once for the part of the space on each side. 

In the case of the isolated points, if the surface integral or line integral round 
this added boundary vanish, we can of course cease to consider these points as 
singular. Suppose Q becomes infinite at the point p = a. Draw a small sphere of 
radius a and also a sphere of unit radius with the point a for centre, and consider 
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the small sphere to be the added boundary. Let dL' be the element of the unit 
sphere cut off by the cone which has a for vertex and the element dL of the small 
sphere for base. Then dL = a 2 dl! and we get for the part of the surface integral 
considered a 2 ffQ d z dL' where Q dIy is the value of Q at the element dL. If then 

Lt T(p _ a)=0 T 2 (p - a)QjzU(p - a) = 

the point may be regarded as not singular. If the limiting expression is finite the 
added surface integral will be finite. If the expression is infinite the added surface 
integral will be generally but not always infinite. Similarly in the case of an 
added line integral if Lt T(fi _ a)=0 T(p - a)QU(p - a) is zero or finite, the added line 
integral will be zero or finite respectively (of course including in the term finite 
a possibility of zero value). If this expression be infinite, the added line integral 
will generally also be infinite. 

This leads to the consideration of potentials which is given in § 9. 

8. Some particular cases of equations (8) and (9) which (except the last) 
have been proved by Tait, are very useful. First put Q = a simple scalar P. Thus 

jPdp = ffVdZVP, (10) 

ffPdZ = JffVPd*. (ID 

If P be the pressure in a fluid -ffPdl, is the force resulting from the pressure on 
any portion and equation (11) shews that -VP is the force per unit volume due to 
the same cause. Next put Qco = Scocr and Vcocr. Thus 

jSdpcr = ffSdUVa, (12) 

jVdpcr = ffV(VdZV . cr) = ffdZSVcr-ffViS dZcr u (13) 

JfSdZcr = fffSVcrds, (14) 

ffVdZcr = fffWo-ds. (15) 

Equations (12) and (14) are well-known theorems, and (13) and (15) will receive 
applications in the following pages. Green's Theorem with Thomson's extension 
of it are, as indeed has been pointed out by Tait particular cases of these equations. 
Equations (14) and (15) applied to an element give the well-known physical 
meanings for S Vcr and VVcr*. The first is obtained by applying (14) to any ele- 
ment, and the second (regarding cr as a velocity) is obtained by applying (15) to 



*[Note added, 1892. Let me disarm criticism by confessing that what follows concerning 
Wcr is nonsense.] 
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the element contained by the following six planes each passing infinitely near to 
the point considered — (1) two planes containing the instantaneous axis of rota- 
tion, (2) two planes at right angles to this axis, and (3) two planes at right angles 
to these four. 

One very frequent application of equation (9) may be put in the following 
form: — Q being any linear function (varying from point to point) of 7?i and V\,R 
being a function of the position of a point 

JDjG(^i,Vi)rfs = -JOjG 1 (/?,v 1 )rfs + JfG(/?, < ffi). (16) 

Potentials. 

9. We proceed at once to the application of these theorems in integration 
to Potentials. Although the results about to be obtained are well-known ones in 
Cartesian Geometry or are easily deduced from such results it is well to give this 
quaternion method if only for the collateral considerations which on account of 
their many applications in what follows it is expedient to place in this preliminary 
section. 

If R is some function of p -p' where p' is the vector coordinate of some point 
under consideration and p the vector coordinate of any point in space, we have 

p VR = - p ,Vfl. 

Now let Q(R) be any function of R, the coordinates of Q being functions of p only. 
Consider the integral JJjQ(R) ds the variable of integration being p (p' being 
a constant so far as the integral is concerned). It does not matter whether the 
integral is a volume, surface or linear one but for conciseness let us take it as a 
volume integral. Thus we have 

^JffQ{R)ds = fff p S7Q(R)ds = -jffF&ROds. 

Now p V operating on the whole integral has no meaning so we may drop the affix 
to the V outside and always understand p' . Under the integral sign however we 
must retain the affix p or p' unless a convention be adopted. It is convenient to 
adopt such a convention and since Q will probably contain some P V but cannot 
possibly contain a p <V we must assume that when V appears without an affix under 
the integral sign the affix p is understood. With this understanding we see that 
when V crosses the integral sign it must be made to change sign and refer only to 
the part we have called R. Thus 

VfffQ(R)ds = -JffViQiROds. (17) 
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Generally speaking R can and will be put as a function of T~ l {p' - p) and 
for this we adopt the single symbol u. Both this symbol and the convention just 
explained will be constantly required in all the applications which follow. 

10. Now let Q be any function of the position of a point. Then the potential 
of the volume distribution of Q, say q, is given by: — 

q = fffuQds, (18) 

the extent of the volume included being supposed given. We may now prove the 
following two important propositions 

V 2 ^ = 4nQ, (19) 

ffSdHV .q = 4nfffQ ds. (20) 

The latter is a corollary of the former as is seen from equation (9) § 6 above. 

Equation (19) may be proved thus. If P be any function of the position of a 
point which is finite but not necessarily continuous for all points 

VfffuPds, 

is always finite and if the volume over which the integral extends is indefinitely 
diminished, so also is the expression now under consideration, and this for the 
point at which is this remnant of volume. This in itself is an important proposi- 
tion. The expression, by equation (17), = - fffv uP ds and both statements are 
obviously true except for the point p'. For this point we have merely to shew 
that the part of the volume integral just given contributed by the volume indef- 
initely near to p' vanishes with this volume. Divide this near volume up into a 
series of elementary cones with p' for vertex. If r is the (small) height and du> 
the solid vertical angle of one of these cones, the part contributed by this cone is 
approximately U(p-p')P p irdu>/3 where P p > is the value of P at the point p'. The 
proposition is now obvious. 
Now since 

V 2 q = V 2 fff u Qds = JffV 2 uQds, 

we see that the only part of the volume integral JJJuQ ds which need be consid- 
ered is that given by the volume in the immediate neighbourhood of p', for at all 
points except p', V 2 w = 0. Consider then our volume and surface integrals only 
to refer to a small sphere with p' for centre and so small that no point is included 
at which Q is discontinuous and therefore VQ infinite. (This last assumes that Q 
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is not discontinuous at p'. In the case when Q is discontinuous at p' no definite 
meaning can be attached to the expression V 2 q.) We now have 

v2 <? = v 2 fff u Q d * = - v HT Vm Q d5 [ b y § 9 ] 

= Vfffu VQ ds - Vffu dL Q, 

equation (9) § 6 being applied and the centre of the sphere being not considered 
as a singular point since the condition of § 7 is satisfied, viz. that LtT(p~ P ')=o T 2 (p- 
p')uU(p - p') = 0. Now putting P above = VQ we see that the first expression, 
viz. V JJJu VQ ds can be neglected and the second gives 

V 2 q = ffVudZQ = 4nQ, 

where Q is the mean value of Q over the surface of the sphere and therefore in 
the limit = Q. Thus equation (19) has been proved. 

When Q has a simple scalar value all the above propositions, and indeed pro- 
cesses, become well-known ones in the theory of gravitational potential. 

We do not propose to go further into the theory of Potentials as the work 
would not have so direct a bearing on what follows as these few considerations. 



Section III. 

*Elastic Solids. 

Brief recapitulation of previous work in this branch. 

11. As far as I am aware the only author who has applied Quaternions to 
Elasticity is Prof. Tait. In the chapter on Kinematics of his treatise on Quater- 
nions, §§ 360-371, he has considered the mathematics of strain with some elab- 
oration and again in the chapter on Physical Applications, §§ 487-491, he has 
done the same with reference to stress and also its expression in terms of the 
displacement at every point of an elastic body. 

In the former he has very successfully considered various useful decompo- 
sitions of strain into pure and rotational parts and so far as strain alone is con- 
sidered, i.e. without reference to what stress brings it about he has left little or 
nothing to be done. In the latter he has worked out the expressions for stress by 
means of certain vector functions at each point, which express the elastic proper- 
ties of the body at that point. 

But as far as I can see his method will not easily adapt itself to the solution of 
problems which have already been considered by other methods, or prepare the 
way for the solution of fresh problems. To put Quaternions in this position is our 
present object. I limit myself to the statical aspect of Elasticity, but I believe that 
Quaternions can be as readily, or nearly so, applied to the Kinetics of the subject. 

For the sake of completeness I shall repeat in my own notation a small part of 
the work that Tait has given. 

Tait shews (§ 370 of his Quaternions) that in any small portion of a strained 
medium the strain is homogeneous and (§ 360) that a homogeneous strain func- 
tion is a linear vector one. He also shews (§ 487) that the stress function is a linear 



*[Note added, 1892. It would be better to head this section "Elastic bodies" since except when 
the strains are assumed small the equations are equally true of solids and fluids. I may say here 
that I have proved in the Proc. R. S. E. 1890-91, pp. 106 et seq., most of the general propositions 
of this section somewhat more neatly though the processes are essentially the same as here.] 
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vector one and he obtains expressions (§§ 487-8) for the force per unit volume 
due to the stress, in terms of the space-variation of the stress. 

Strain, Stress-force, Stress-couple. 

12. This last however I give in my own notation. His expression in § 370 
for the strain function I shall throughout denote by x- Thus 

XO) = to - S to V . r\, (1) 

where n is the displacement that gives rise to the strain. 

Let^ consist of a pure strain if/ followed by a rotation q{ )q~ l as explained in 
Tait's Quaternions, § 365 where he obtains both q and if/ in terms of X- Thus 

XO) = q if/co q~ l . (2) 

When the strain is small if/ takes the convenient form x where x stands for the 
pure part of x so that 

XO) - to - \ScoV . rj - ^ViScor/i, (3) 

by equation (3) §2 above. Similarly q( )q~ l becomes V6{ ) where 29 = Wn. 
The truth of these statements is seen by putting in equation (2) for q, 1 + 6/2 and 
therefore for q~ x , 1 - 6/2 for if/, x and then neglecting all small quantities of an 
order higher than the first. 

13. Next let us find the force and couple per unit volume due to a stress 
which varies from point to point. Let the stress function be (p. Then the force on 
any part of the body, due to stress, is 

ff(f>dZ = fff<f>Ads, 

by equation (9) § 6. Thus the force per unit volume = <pA, for the volume consid- 
ered in the equation may be taken as the element ds. 
Again the moment round any arbitrary origin is 

ffVpctdZ = JffVpKtKx ds + fffVp^V, ds, 

by equation (9) § 6. The second term on the right is that due to the force cf>A just 
considered, and the first shews that in addition to this there is a couple per unit 
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volume = V£(p£, = twice the "rotation" vector of (p. Let then tp be the pure part 
and Ve( ) the rotatory part of tp. Thus 
Force per unit volume 



(f>A = (/)A- We, (4) 



Couple per unit volume 



= VCK = 2e. (5) 

These results are of course equivalent to those obtained by Tait, Quaternions, 
§§ 487-8. 

The meanings just given to r\, x, X, iff, q, tp, tp and e will be retained throughout 
this Section. In all cases of small strain as we have seen we may use iff or x indif- 
ferently and whenever we wish to indicate that we are considering the physical 
phenomenon of pure strain we shall use if/, x being regarded merely as a function 
of x- We shall soon introduce a function m which will stand towards (p somewhat 
as if/ towards x and such that when the strain is small m = (p. 

It is to be observed that (poo is the force exerted on a vector area, which when 
strained is to, not the stress on an area which before strain is to. Similarly in 
equations (4) and (5) the independent variable of differentiation is p + n so that 
strictly speaking in (4) we should put (p p+T] A - V p+t] Ve. In the case of small strain 
these distinctions need not be made. 

Stress in terms of strain. 

14. To express stress in terms of strain we assume any displacement and 
consequent strain at every point of the body and then give to every point a small 
additional displacement Si] and find in terms of if/ and <p the increment JJJSw dso 
in the potential energy of the body, w d% being the potential energy of any ele- 
ment of the body whose volume before strain was d%. Thus 

JJJSw dso = (work done by stresses on surface of portion considered) 

- (work done by stresses throughout volume of same portion). 

Thus, observing that by § 12 the rotation due to the small displacement 6n is 
V p+ri VSri/2, we have 

fffSw d% = -ffS 5r]cpd'L + fffS S^tp^Vi ds + fffSe p+n VSrids. 

The first of the terms on the right is the work done on the surface of the portion of 
the body considered; the second is -(work done by stress-forces P+;? A); and the 
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third is -(work done by stress-couples 2e). Thus converting the surface integral 
into a volume integral by equation (9) § 6 above 

fffSwd% = -JffS 6n x cf) p+v Vi ds + fffs 6rji Ve^Vi ds 
= -fffS6 m p+r] V l ds. 

Limiting the portion of the body considered to the element ds> we get 

Sw = -mS Srii ~4> p+ ^\ (6) 

where m is put for ds/ds and therefore may be put by § 3a and § 3 above in the 
various forms 

6m = S & && Sx€iX&X& ( 6a ) 

= S{i &6 Sx'hx'bx'b (fib) 

= Ste 2 bSMMi»l'& (6c) 

= 5V 1 V 2 V 3 5(p + i7) 1 (p + i7)2(p + i7)3. (6d) 

It is to be observed that since the rotation-vector e of (p does not occur in 
equation (6), e and therefore the stress-couple are quite arbitrary so far as the 
strain and potential energy are concerned. I do not know whether this has been 
pointed out before. Of course other data in the problem give the stress-couple. In 
fact it can be easily shewn that in all cases whether there be equilibrium or not 
the external couple per unit volume balances the stress-couple. [Otherwise the 
angular acceleration of the element would be infinite.] Thus if M be the external 
couple per unit volume of the unstrained solid we have always 

M + 2me = 0. (7) 

In the particular case of equilibrium F being the external force per unit volume of 
the unstrained solid we have 

F + m(0 p+n A - V p+J1 Ve) = 0. (8) 

The mathematical problem is then the same as if for ¥/m we substituted ¥/m + 
V p+rj V(M/2m) and for M, zero. In the case of small strains m may be put = 1. 
In this case then the mathematical problem is the same as if for F we substituted 
F + VVM/2 and for M zero. 

Returning to equation (6) observe that 

6xo) = -Sxoj p+v V . 8r\ 



28 ELASTIC SOLIDS. [§ 14. 

[which is established just as is the equation xco = -ScoV . (p + rj)], so that chang- 
ing co which is any vector into;^ 1 ^ 

Sx.X~ l u = Sco p+r] V .6q. 
Therefore by equation (6) of this section and equation (4) § 3 above we have 

Sw = -mS5xx~ l £H- 

We must express the differential on the right of this equation in terms of dip 
and 8q, the latter however disappearing as we should expect. Now^w = qif/coq~ l 
[equation (2) § 12] so that remembering that 6 . q' 1 = -q' 1 5qq~ l 

5xco = Sq if/coq~ l - qif/coq' 1 Sqq~ l + q8\f/ coq~ l 
= 2VV 5qq~ l -x^o + qSif/coq' 1 , 
5w/m = -2S . V6qq~ 1 . £#£ - Sq6if/x~H <f l ~4>^ 

or since V £cf>£ = 0, tb being self-conjugate 

8w = -mS Sif/x l £><l ty^l- 

This can be put into a more convenient shape for our present purpose. First put 
for x' 1 its value \]/~ l q~ l { )q and then apply equation (6a) § 3 above, putting for 
the 4> of that equation a _1 ( )q and therefore for the <p' ', q( )q~ 1 . Thus 

8w = -mSSifftlf- l £i7jt, (9) 

where mco = q~ l cb{qcoq~ l )q. (10) 

The physical meaning of this last equation can easily be shewn. Suppose 
when there is no rotation that cf> = m' '. Then it is natural* to assume — in fact it 
seems almost axiomatic — that the superimposed rotation q( )q~ x should merely 
so to speak rotate the stress along with it. Thus if co is some vector area before 
the rotation which becomes co' by means of the rotation 

qvj'co . q~ l = cbco'. 

But co' = qcoq~ l , so that 

tu'co = q' i cf>(qcoq' i )q = vjco. 



'Observe that we do not make this assumption. We really shew that it is true. 
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Thus we see that m is as it were cp with the rotation undone. 

Before proceeding further with the calculation let us see what we have as- 
sumed and what equation (9) teaches us. The one thing we have assumed is that 
the potential energy of the body can be taken as the sum of the potential energies 
of its elements, in other words that no part of the potential energy depends con- 
jointly on the strains at P and Q where P and Q are points separated by a finite 
distance. This we must take as an axiom. By it we are led to the expression for 6w 
in equation (9). This only involves the variation of the pure strain if/ but not the 
space differential coefficients of if/. This is not an obvious result as far as I can see 
but it is I believe always assumed without proof. 

We may now regard w as a function of iff only. Therefore by equation (7) § 5 
above 

6w = -S dif/t; ^qw g, 

therefore by equation (9) of this section 

S Sif/ £ ^aw £ = mS Sif/if/' 1 r,wC ■ 

In equation (6a) § 3 above putting cp = if/' 1 the right-hand member of this equa- 
tion becomes mS 5if/ £ mif/~ l £. Now putting in equation (6a) §30 = mif/' 1 this 
member becomes mS^dif/if/' 1 ml, or mS 6if/£if/~ l m£. Thus we have 

2S 5if/^aw£ = mS Sif/ ^{mif/' 1 + iff^m)^. 

Now xjjif/~ l + ip~ x vj is self-conjugate. Hence by § 4 above 

m{vjif/~ x + iff' 1 m) = 2 iA qw. (11) 

This equation can be looked upon as giving m in terms of the strain. We 
can obtain m however explicitly for i//~ 1 tj is the conjugate of miff' 1 . Hence from 
equation (11) (because (if/^m + mif/~ l )/2 is the pure part of mif/~ l ) 

mmif/'^d) = ^(jwd) + V6co 

where is a vector to be found. Changing to into if/to 

m mio = ^Qw iff co + V6 iff to. 



*By putting 8tff = a>S ( )a>' + co'S ( )a> in this equation, equation (11) can be deduced but as this 
method has already been applied in § 4 I give the one in the text to shew the variety of Quaternion 
methods. [Note added, 1892. If we use the theorem in the foot-note of § 4, equation (11) follows 
at once.l 
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Now vj being self-conjugate V^mZ, = 0. Hence 

v^awH = -vzvew = esm - Msec 



whence 
or e = (ilf + S^ilfC\r X VH*QwC 



by equation (6a) § 3 above. Thus finally 
mmco = ^qw if/to + VOif/to 



(12) 



or mmco = ^Qwif/co - Vif/co{\J/ + S£\\J/£,\) Vif/^Qw g. 



(13) 



This completely solves the problem of expressing stress in terms of strain in 
the most general case. 

15. into + Ve'io, where, as we saw in last section, e' is perfectly arbitrary so 
far as the strain is concerned, is the force on the strained area to due to the pure 
strain iff. And again 

tbto + Ve"to or q{mq~ l wq)q~ x + Ve"a> 

is that due to the strain qif/( )q~ { or x- 

To find the force on an area which before strain was to let its strained value 
after iff has taken place be to. Then by equation (4), § 145 of Tait's Quaternions, 
ma>o = ip to. Hence 

Required force = into + Ve'co = ,/,Qwto + V6co + mVe , i//~ 1 a>o 

by equation (13) of last section. If the rotation now take place this force rotates 
with it so that the force on the area which was originally to Q is after the strain x 
or qiff{ )q- 1 

tcoq = q^Qw to Q q~ l + qV0to Q . q~ l + V(eqiff~ coq . q~ l ) (14) 

where e = mqe'q 1 and 

e=2 tne stress-couple per un. vol. of unstrained body. 
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This force then is a linear vector function of a>o, but in general even when 
e = it is not self-conjugate. When both e and the rotation are zero we see that 
the rotation vector of t is 9 given by equation (12) of last section. 

The stress-force can be shewn as in § 13 to be tA per unit volume of the 
unstrained body. Thus since the corresponding stress-couple is 2e the moment 
exerted by the stresses on any portion of the body round an arbitrary origin is 

fffV(p + jflTiVj rfs + 2fffeds . 

But this moment may also be put in the form 

ffV(p + T])TdY , 

or JffV^dso + fffV(p + iflTxVj rfs + fffVij^i ds , 

by equation (9) § 6 above. Comparing these results 

2e = V{T{ + VT h TV l \ 

or , > (14a) 

2e=Vp' 1 TV 1 , J 

in the notation of next section. This equation may also be deduced from equa- 
tion (15/) below but not so naturally as above. 

The equations of equilibrium. 

16. We require equation (9) § 14 above to prove the statement that no space- 
variations of if/ or q are involved in w. It is also required to shew that the fact 
that q also is not involved in w is a mathematical sequence of the assumption that 
the potential energy of a solid is the sum of the potential energies of its elements. 
Assuming these facts however we can arrive at the equation of stress (1 1) § 14 in a 
different way from the above. We shall also obtain quite different expressions for 
m, t, &c. and most important of all we shall obtain the equations of equilibrium 
by obtaining r explicitly in terms of the displacement and its space derivatives. 
From § 12 above we have 

X<o = qif/ioq' 1 , x'w = *p{q~ l toq), 

X 'X = V = Y, (15) 

(say) as in Tait's Quaternions, § 365. Thus ^co = V \S coV ' 2 S p\p' 2 where p' is 
put as it will be throughout this section for p + q, the vector coordinate after 
displacement of the point p. From this as we have seen in the Introduction we 
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deduce that the coordinates of *P are the A, B, C, a, b, c of Thomson and Tait's 
Nat. Phil. App. C. 

We may as do those authors regard w as a function of X P. Thus: — 

= S dx'x^Gw £ - Sx' Sx &QW £. 
By equation (4a) § 3 above, each of the terms in this last expression 

= -S SrjiXvQwVu 

6w = -2S SrjiXyQwVi, (15a) 

Comparing this with equation (6) § 14 and putting in both 5rji = co'Scop' where 
co', co are arbitrary constant vectors, we get 

mSco' 0to = ISto'xvQwx'w- 

Hence, since co' is quite arbitrary 

imp = 2xy(jwx'- (15a) 

From equation (10) § 14 above which defines m we see that 

mm = lif/xyQwiff, (15c) 

which we should expect since we have already seen that m is the value of cp when 
there is no rotation and therefore x = x' = &• Now since 

S Sifr^aC = S &¥fra{ = S(6i/fiff + ijr6ifr)frag, 

we deduce by any one of the processes already exemplified that 

^a = wQi/f + tffvQ, 

where of course the differentiations of >pQ must not refer to if/. We see then from 
equation (15c) that 

m{m^ x + i/r'tzr) = 2^qw, 

which is equation (11) § 14. 

Our present purpose however is to find the equations of equilibrium. Let coq 
be the vector area which by the strain ;f becomes co. Thus* as in last section 

mco Q = x' 1 ^- (15c) 



*See § 83 below. 
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Further let 2e be the stress-couple per unit volume of the unstrained solid so 
that 2e/m is the same of the strained solid. As we know, e is quite independent 
of the strain. By § 13 we see that the force tco on the area which before strain 
was coq is (pco + Veco/m. Therefore 

tcoq = 2xwawco Q + Vey~W (15/) 

We saw in the last section that the force per unit volume of the unstrained solid 
is tA and the couple 2e. Hence 

F + tA = 0, (I5g) 

M + 2e = 0, (15/0 

are the equations of equilibrium, where F and M are the external force and couple 
per unit volume of the unstrained solid. All that remains to be done then is to 
express t in terms of e and the displacement. Putting 

p + 7]=p', (15/) 

as already mentioned, we have 

Xoj = -p\ScoV u (15;) 

so that by equation (6/0 § 3a above we have 

x'-'oj = -3Vp' l p' 2 S(oV l V 2 /SV 1 V 2 V i Sp\p' 2 p' 3 . (15k) 

Therefore by equation (15/) 

to) = -2p\SV lw awoj- 3VeVp\p' 2 S(xN l V 2 ISV l V 2 V^Sp\p' 2 p' y (15/) 

It is unnecessary to write down what equation (I5g) becomes when we substitute 
for t, changing e into -M/2 by equation (15/0. In the important case however 
when M = 0, the equation is quite simple, viz. 

F = 2p\SV lw QwA. (15m) 

Addition to § 16, Dec, 1887 (sent in with the Essay). [The following consid- 
erations occurred just before I was obliged to send the essay in, so that though I 
thought them worth giving I had not time to incorporate them in the text. 

It is interesting to consider the case of an isotropic body. Here w is a function 
of the three principal elongations only and therefore we may in accordance with 
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§ 14 and § 15 suppose it a function of a, b, c or in accordance with § 16 of A, B, 
C where 

a = -S£if/£, b = -S^ 2 i;/2, c = -S£if/ 3 £/3, (A) 

A = -SfPi;, B = -Sp¥ 2 £/2, C = -S<pF 3 £/3. (B) 

Let us use x, y, z, X, Y, Z for the differential coefficients of w with respect to 
a, b, c, A, B, C respectively. Thus 

dw = xda + ydb + zdc = -S difr £(x + yif/ + z^ 1 )^ 

as can be easily proved by means of equation (6a) § 3 above. But 

dw = -S dif/^Qw g. 

^aw = x + yif/ + zfi 2 , (C) 

by § 4 above. Similarly we have 

v(jw = X + Y X ¥ + Z*¥ 2 . (D) 

(Notice in passing that to pass to small strains is quite easy for ^qw is linear and 
homogeneous in if/ - 1 so that ^qw = x + yip where y is a constant and x + y a 
multiple of a - 3.)* 

From equation (C) we see that in equation (12) § 14 is zero and therefore 
that from equation (13) 

mm = p(jw if/ = xifj + yif/ 2 + ziff 3 . (E) 

Again from equation (14) § 15 

to) = xqa>q~ l + yxoj + zxfto) + Vex'~ l oj- (F) 

Similarly from equations (15&) and (15/) § 16 we may prove by equation (D) that 

\m~4> = X XX ' + YCf/) 2 + Z(xx')\ (G) 

and r = 2X X + 2Y XX ' X + ^XX' XX' X + Vex"\ )■ (#) 

If we wish to neglect all small quantities above a certain order the present 
equations pave the way for suitably treating the subject. I do not however pro- 
pose to consider the problem here as I have not considered it sufficiently to do it 
justice.] 



"[Note added, 1892. In the original essay there was a slip here which I have corrected. It 
was caused by assuming that if/ instead of ifi - 1 was small for small strains. In the original I said 
"where y is a constant and x is a multiple of a."] 
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Variation of Temperature. 

16a. The w which appears in the above sections is the same as the w which 
occurs in Tait's Thermo-dynamics, § 209, and therefore all the above work is true 
whether the solid experience change of temperature or not. w will be a function 
then of the temperature as well as if/. To express the complete mathematical 
problem of the physical behaviour of a solid we ought of course instead of the 
above equations of equilibrium, to have corresponding equations of motion, viz. 
equations (15/z), (15/) and (instead of (15g)) 

Dp' = F + rA, (15n) 

where D is the original density of the solid at the point considered. Further we 
ought to put down the equations of conduction of heat and lastly equations (16/) 
and (166?) below. 

We do not propose to consider the conduction of heat, but it will be well to 
shew how the thermo-dynamics of the present question are treated by Quater- 
nions. 

Let t be the temperature of the element which was originally d% and E ds its 
intrinsic energy. Let 

6H = -S 5if/£M£ + NSt, 

where M, N are a linear self-conjugate function and a scalar respectively, both 
functions of if/ and t. Here 6H d% is the heat required to be put into the element 
to raise its temperature by 5t and its pure strain by 5if/. Now when t is constant 
5H must be a perfect differential so that we may put 

M = t+Qf, 

where / is some function of if/ and t. Thus 

6H = -tS dif/^Qft + NSt. (16) 

Now we have seen in § 14 that the work done on the element during the incre- 
ment 5ifr, divided by d% 

= -mS 6if/if/~ l £tzt£. 

Thus by the first law of Thermo-dynamics 

8E = JSH-mS Si// if/' 1 t^w^ 

= JNSt- JtS Sif/^aft ~ mS Sif/if/~ l £tj£, 
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where / is Joule's mechanical equivalent. Thus 

JN = dE/dt, (16a) 

and \m(m\]j~ l + if/~ l m) = ^Qw, (16b) 

where w = E - Jtf. (16c) 

To apply the second law we go through exactly the same cycle as does Tait in 
his Thermo-dynamics, § 209, viz. 

(if/, t)(if/ + 6if/, t)(if/ + 6if/,t + 6t)(if/, t + 8i)(if/, i). 

We thus get* 

j (dm _j _ x dvj\ d,/,(jw dw 



h^S = - 1 m\-r +r -\ = -— = -, a - 



or // = -dw/dt (I6d) 

the arbitrary function of t being neglected as not affecting any physical phe- 
nomenon. Substituting for w from equation (16c), 

Jtdf/dt = dE/dt = JN (16e) 

by equation (16a). Thus from equation (16) 

SH = tSf, (16/) 

so that in elastic solids as in gases we have a convenient function which is called 
the "entropy". Thus the intrinsic energy E d%, the entropy f d% and the stress m 
have all been determined in terms of one function w of if/ and t which function is 
therefore in this general mathematical theory supposed to be known. 

If instead of regarding w (which with the generalised meaning it now bears 
may still conveniently be called the "potential energy" per unit volume) as the 
fundamental function of the substance we regard the intrinsic energy or the en- 
tropy as such it will be seen that one other function of if/ must also be known. For 
suppose / the entropy be regarded as known. Then since dw/dt = -Jf 

w = W-Jffdt, (I6g) 



*[By assuming from the second law that the work done by the element in the cycle, i.e. the 
sum of the works done by it during the first and third steps is J8t ft multiplied by the heat absorbed 
by the element in the third step. Note added, 1893.] 
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where the integral is any particular one and W is a function of if/ only, supposed 
known. Again 

E = w + Jtf 
or E = W + J(tf- Jfdt). (16/0 

Thus all the functions are given in terms of / and W. Similarly if E be taken as 
the fundamental function 

w = t(W - jEdtft 2 ), (16/) 

Jf = Eft + jEdtft 2 -W, (16/) 

where as before the integral is some particular one and W is a function of if/ only. 

Small strains. 

17. We now make the usual assumption that the strains are so small that 
their coordinates can be neglected in comparison with ordinary quantities such as 
the coefficients of the linear vector function ^,qw. We can deduce this case from 
the above more general results. 

To the order considered q = 1 so that by equation (10) § 14 above 

<f> = VJ. 

We shall use the symbol w rather than (p for the same reasons explained in § 13 
above as induce us to use if/ rather than^. 

Remembering that aj and a> are now the same and that if/ may be put = 1 so 
that V^Qwif/^ = V£^Qwg = and therefore of equation (12) = 0, we have 
from equation (13), 

vju) = ^(jwco. (17) 

Of course we do not require to go through the somewhat complicated process 
of § 14, § 15 to arrive at this result. In fact in equation (6) § 14, we may put m = 1 
and P+?7 V = V so that 

5w = -S dr/i 0Vi = -S dr/i mV\, 

and therefore by § 3 above 

8w = -S Sx£ut£ = -S 5if/ £,vj£,. 

But by § 5 above 

5w = -S 6if/ ^qw £, 

and therefore by § 4 we get equation (17). 
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18. It is convenient here to slightly change the notation. For if/, x we shall 
now substitute if/ + 1, y + 1 respectively. This leads to no confusion as will be 
seen. 

With this notation the strain being small the stress is linear in if/ i.e. ^qw is 
linear and therefore w quadratic. Now for any such quadratic function 

w = -Sif/^Qwa2, (18) 

for we have by §5, 

dw = -S dif/ £#Qw £,. 
Put now if/ = nif/'. Then because ^qw is linear in if/ 

where ^(j'w is put for the value of ^qw when the coordinates of if/' are substituted 
for those of if/. Thus keeping if/' constant and varying n, 

dw = -ndnS if/' ^(j'w g, 

whence integrating from n = to n = 1 and changing if/' into if/ we get equa- 
tion (18). 

From equation (18) we see that for small strains we have 

w = -Sif/^/2. (19) 

Now w is quadratic in if/ and therefore also quadratic in m, so that regarding w 
as a function of m we have as in equation (18) 

w = -Sm^(jw£/2, 

so that by equation (19) and § 4 above 

if/ = m aw. (20) 

All these results for small strains are well-known in their Cartesian form, but 
it cannot be bias that makes these quaternion proofs appear so much more natural 
and therefore more simple and beautiful than the ordinary ones. 
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19. Let us now consider (as in § 16 is really done) w as a function of the 
displacement. Now w is quadratic in if/, and iff is linear and symmetrical in Vi 
and rji. In fact from equation (3) § 12 above, remembering that the x of that 
equation is our present if/ + 1 , we have 

2iffto = -r\\ ScoVi - Vi Scorji. 

Therefore we may put 

w = w(77i,Vi,7/ 2 ,V 2 ), (21) 

where w(a,j3, y, 6) is linear in each of its constituents, is symmetrical in a and ft, 
and again in y and 5, and is also such that the pair a, fi and the pair y, 6 can 
be interchanged. [This last statement can be made true if not so at first, by sub- 
stituting for w{a,/5,y, 8), w(a,fi,y,5)/2 + w(y,S,a,{3)/2 as this does not affect 
equation (21).] Such a function can be proved to involve 21 independent scalars, 
which is the number also required to determine an arbitrary quadratic function 
of if/, since iff involves six scalars. 

Thus we have the two following expressions for Sw, which we equate 

-S Sifj^mt; = w(6r] U Vi,t]2, V 2 ) + w(tji,V u Stj 2 , V 2 ) 
= 2w(Sr} 1 ,V l ,rj2,V 2 ), 

or* by § 3 above, 

-S SijiwVi =2w(&7i,Vi,77 2 ,V 2 ). 

Now let us put 6rj = co'Scop where co' and to are arbitrary constant vectors. We 
thus get 

Sco' lino = -2w((o',a>, rji, Vi) = 2Sto' £w(£, co, r\\, Vi). 

Whence since to' is quite arbitrary, 

etg> = 2£w(£,o>,77i,Vi). (22) 

The statical problem can now be easily expressed. As we saw in § 14, equa- 
tions (7) and (8), it is simply 

F + VVM/2 + gtA = 0, (22a) 



*[Note added, 1892. Better thus: — by § 3 above, 

-S 6if,^; = 2w{6iff^;, i;i,Vi) = -25 fi ^fwtfi.f.iji.Vi) 

therefore by § 4, tltoj = 2(w(£, co, 771, Vi) 

for (w(£, to, 771, Vi) regarded as a function of co is clearly self conjugate.] 
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throughout the mass; and at the surface 

F s - VU rfZ M/2 -vrUdL = 0, (22b) 

where F, M are the given external force and couple per unit volume and ¥ s is 
the given external surface traction per unit surface. Substituting for vj from equa- 
tion (22) 

F + VVM/2 + 2£w(£,A,i/i,Vi) = 



(23) 
F 5 - VU d!M/2 -2£w(£,U d^TjuVi) = 0.} 

Isotropic Bodies. 

20. The simplest way to treat these bodies is to consider the (linear) relations 
between tu and iff. 

In the first place notice that iff can always be decomposed into three real elon- 
gations (contractions being of course considered as negative elongations). Thus i 
being the unit vector in the direction of such an elongation, 

iff to = -Lei Si to. 

The elongation -eiSico will cause a stress symmetrical about the vector i, i.e. a 
tension Ae in the direction of i and a pressure Be in all directions at right angles; 
A and B being constants (on account of the linear relation between m and iff) 
independent of the direction of i (on account of the isotropy of the solid). This 
stress may otherwise be described as a tension (A + B)e in the direction of i and a 
hydrostatic pressure Be. Thus 

vjco = -(A + S)E ei Si to - Bco ~Le 
= (A + B)if/co + BcoS^iff^. 

To obtain the values of A and B in terms of Thomson and Tait's coefficients k 
and n of cubical expansion and rigidity respectively; first put 

if/to = eu> and mo) = 3keto, 

and then put 

iff a> = VAtopi and mco = 2nVAcop., 
A and ii being any two vectors perpendicular to each other. We thus get 

A + B = 2n, B = -(k - 2n/3), 
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whence 

wco = 2n i//cl> - (k - 2n/3)a>S £i//£. (24) 

From this we have 

if/to = mco/2n + coS £i[/£(k - 2n/3)/2n, 
but from the same equation 

S^ = 3kS£M; 

1 Wl-JLV 



il/io = — tuoj + coS Cxnt. (25) 

v 2n \6n 9k) s s 

Equation (24) gives stress in terms of strain and (25) the converse. 



21. We can now give the various useful forms of w for isotropic bodies for 
from equation (19) § 18, 

w = -Sif/ gTcrg/2. 

Therefore from equations (24) and (25) respectively 

w = -n{Uf + i(k - 2n/3)S £M\ S bWi, (26) 

w = ~(^-i(^-^^i5^ 6 . (27) 

Therefore again from § 3 above and from equation (26), 

w = -nSTju/rVx + Uk - 2n/3)(SVi]) 2 



or 



2w(t7i,Vi,^ 2 ,V 2 ) = nSViT] 2 SV 2 r]i + nSViV 2 S 771/72 

+(k-2n/3)SV 1 T ll SV 2 ri2- 



(28) 



Hence from equation (22) 

vtoj = -nSa>V . t] - nViSrjito - (m - n)coSVrj, (29) 

where m is put for k + n/3. This last could have been deduced at once from 
equation (24) by substituting for if/a>. 
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Thus the equations (23) for the statical problem are 

F + VVM/2 = nV 2 ri + mVS V17, (30) 

-F s + VUdZM/2 = nSUdZV .rj + nViS^UdX (31) 

+ (m-n)Ud'LSVr]. 

We now proceed to apply these results for small strains in isotropic bodies 
to particular cases. These particular cases have all been worked out by the aid of 
Cartesian Geometry and they are given to illustrate the truth of the assertion made 
in the Introduction that the consideration of general problems is made simpler by 
the use of Quaternions instead of the ordinary methods. 

Particular integral of equation (30)*. 

22. Since from equation (30) (F being put for simplicity instead of F + 
VVM/2) we have 

nV 2 r] = F -mVSVr] 

we obtain as a particular case by equations (18) and (19) § 10, 

Annr] = JJJu(F - mVSVrj) ds, 

where u has the meaning explained in § 9, and the volume integral extends over 
any portion (say the whole) of the body we may choose to consider. To express 
JJJuVSVrjds as a function of F put in this term u = -\Up\V\ where p is taken 
for the p - p' of § 9, and apply equation (9) § 6. Thus 

Annri = fffuFds + ^mfffl/p^iVSVrjds 

= fffuFds - ±mfffupV 2 SVrids + a surf. int. 

= fffuFds fffUpS VFds + the surf. int. 

JJJ 2(m + n) JJJ 

for by equation (30) S VF = (m+n)V 2 S Vrj. Now (in order to get rid of any infinite 
terms due to any discontinuity in F) apply equation (9) § 6 to the second volume 
integral. Thus 

4nnr] = fffuF ds + fffSFV .Upds+ a surf. int. 

JJJ 2(m + n) JJJ 



"[Note added, 1892. For a neater quaternion treatment of this problem see Phil. Mag. 
June, 1892, p. 493.] 
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The surface integral may be neglected as we may thus verify. Call the volume 
integral Annrj' . Thus 

AnnV 2 ri' = 4nF + (jfjfSFV . V 2 Upds 

2(m + n) JJJ 

2 

AnmVSVrj' = - fffSFV . Vwrfs + (jfjTSFV . VSVUpds, 

n JJJ 2n(m + n) JJJ 

so that putting Vf/p = -2w we get 

nVV + mVSVr]' = F, 
whence we have as a particular solution of equation (30) rj = rj' or 

r,= fffuF ds + fffS FV . Up ds. (32) 

' 4nn JJJ 8nn(m + n) JJJ H 

This is generally regarded as a solution of the statical problem for an infinite 
isotropic body. In this case some law of convergence must apply to F to make 
these integrals convergent. Thomson and Tait {Nat. Phil. § 730) say that this 
law is that Fr converges to zero at infinity. This I think can be disproved by an 
example. Put* Fr = r~ a A where A is a constant vector and a a positive constant 
less than unity. Equation (32) then gives for the displacement at the origin due to 
the part of the integral extending throughout a sphere whose centre is the origin 

and radius R 

m + 3n R l a 

rj = A. 

3n(m + n) 1 - a 

Putting R = oo, rj also becomes oo. The real law of convergence does not seem 
to me to be worth seeking as the practical utility of equation (32) is owing to the 
fact that it is a particular integral. 

The present solution of the problem has only to be compared with the one 
in Thomson and Tait's Nat. Phil. §§ 730-1 to see the immense advantage to be 
derived from Quaternions. 

It is easy to put our result in the form given by them. We have merely to 
express S F V . Up in terms of F and r 2 S F V . Vu where r is put for the reciprocal 
of u. Noting that 

Vw = -u 3 p, Up = up, SFV . p = -F, 



"[Note added, 1892. This is not legitimate since it makes F = oo for r = 0. The reasoning 
is rectified in the Phil. Mag. paper just referred to by putting F = from r = to r = b and 
Fr = r~"A from r = b to r - oo.] 
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we have at once 

SFV . Up = -u¥-u 3 pSFp 
SFV .Vu = u 3 F + 3u 5 pSFp 

therefore eliminating pSFp, 

SFV . Up = -uF2/3 - r 2 SFV . V«/3, 
7] = {24nn(m + n)r l fffds {2(2m + 3n)uF - mr 2 SFV . Vw}, (33) 

which is the required form. 

23. Calling the particular solution rf as before and putting 

rj = rj' + t]" 

the statical problem is reduced to finding rj" to satisfy 

nVV + mVSVr]" = 

and the surface equation either 

j]' + j]" = given value, 
i.e. 77" = given value, 

or tjjUcTL = given surface traction, 

i.e. by equation (29) § 21 above, 

nS U dZV . r]" + nViSr]" U dZ + (m - n)U d~LS Vr]" = known value. 

This general problem for the spherical shell, the only case hitherto solved, I 
do not propose to work out by Quaternions, as the methods adopted are the same 
as those used by Thomson and Tait in the same problem. But though each step of 
the Cartesian proof would be represented in the Quaternion, the saving in mental 
labour which is effected by using the peculiarly happy notation of Quaternions 
can only be appreciated by him who has worked the whole problem in both no- 
tations. The only remark necessary to make is that we may just as easily use 
vector, surface or solid, harmonics or indeed quaternion harmonics as ordinary 
scalar harmonics. 
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Orthogonal coordinates. 

24. It is usual to find what equations (22a) of § 19 and (3) of § 12 become 
when expressed in terms of any orthogonal coordinates. This can be done much 
more easily by Quaternions than Cartesian Geometry. Compare the following 
investigation with the corresponding one in Ibbetson's Math. Theory of Elasticity, 
Chap. V. 

Let x, y, z be any orthogonal coordinates, i.e. let x = const., y = const., 
z = const., represent three families of surfaces cutting everywhere at right angles. 
Particular cases are of course the ordinary Cartesian coordinates, the spherical 
coordinates r, 6, (p and the cylindrical coordinates r, (p, z. Let D x , D y , D z stand 
for differentiations per unit length perpendicular to the three coordinate surfaces 
and let A, p, v be the unit vectors in the corresponding directions. Thus 

V = AD X + nD y + vD z . 

Thus, using the same system of suffixes for the D's as was explained in connection 
with V in § 1 , 

(pt± = D xl <p l A + D yl <p l n + D zl <p l v, (34) 

or 

0A = D x {(pA) + D y (<pfi) + D z {(pv) - (p{D x A + D y p + D z v). (35) 

25. Now to put equation (22a) § 19 into the present coordinates all that is 
required is to express crA in terms of those coordinates. Let the coordinates of m 
be PQRSTU. Thus from equation (35) we have 

mA = D X (PA + U/X + Tv) + D y (UA + Qn + Sv) + D Z (TA + S/J+ Rv) 
- m{D x A + D y fi + D z v). 

The first thing then is to find D X A, D x /u, D y A &c. 

Let p 2 , P3 be the principal curvatures normal to x = const., i.e. (by a well- 
known property of orthogonal surfaces) the curvatures along the lines of inter- 
section of x = const., with z = const., and y = const. p 2 , P3 will be considered 
positive when* the positive value of dx is on the convex side of the corresponding 
curvatures. Similarly for qj,q\r\r<i. Thus for the coordinates r, 0, <p; 

P2 = P3 = l/r, q 3 = cot6/r, q l =r 1 = r 2 = 0. 



"[Note added, 1892. This is contrary to the usual convention. 
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Again for r, <p, z', pi = 1 If and the rest are each zero. 
With these definitions we see geometrically that 

D X A = -/xqi-vru D x fx = Aq u D x v = Ar x . (36) 

Similarly for D y A, D y /u, D Z A &c. Thus 

gtA = A(D X P + D y U + D Z T) + //( ) + v( ) 

+ ^(-2^2 - fl/>3 + 7>i + U qi ) +//() + v( ) 

+ ^{(P2 + Ps)^ + (?3 + ^i)f/ + fa + r 2 )r} + ju{ } + v{ }, 

or 

wA = A{D X P + D y U + DT + P(p 2 + p 3 ) - Qp 2 - Rp 3 

+ r(2n + r 2 ) + £% 3 + 2 ?1 )} +//{} + v{ }. (37) 

26. The other chiefly useful thing in transformation of coordinates in the 
present subject is the expression for the strain function if/ in terms of the coor- 
dinates of displacement. Let u, v, w be these coordinates. Now by equation (3) 
§ 12, remembering (§ 18) that if/ = x - 1 we have 

-2(/^a» = ScoV . rj + ViScorji, 

whence -2if/A = AS AD x rj + /jS AD y r] + vS AD z r] - D x rj. 

But D x rj = AD x u + piD x v + vD x w + uD x A + vD x fi + wD x v 

= A(D x u + vq\ + wti) + p.{D x v - uq\) + v(D x w - ur{). 

Similarly 

D y T] = A(D y u - vpi) + n{D y v + w r 2 + wp 2 ) + v(Z> y w - vr 2 ) 
Z) z ?7 = A(D z u - wp 3 ) + fi(D z y - wg 3 ) + v(D z w + up 3 + vg 3 ) 
2t/>vt = 2A(D x u + v^i + wri) + pi(D y u + D^v - w^i - v/? 2 ) 

+ v(D z u + D x w - ur\ - wpi). (38) 

But with Thomson and Tait's notation for pure small strain 

2if/A = 2Ae + /jc + vb, 
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*e = D x u + vq\ + wr\ 

f = D y v + wr 2 + up 2 

g = D z w + up 3 + vq 3 

a = D y w + D z v - wq 3 - vr 2 

b = D z u + D x w - ur\ - wp 3 

c = D x v + D y u - vp 2 - uq\ . 



(39) 



Thus we have efgabc in terms of the displacement and we have already 
in equation (24) § 20, which expresses vj in terms of if/, found the values of 
PQRSTU in terms of e, &c. Finally the expression in equation (37) for wA 
gives us the equations of equilibrium in terms of P, &c. Thus we have all the 
materials for considering any problem with the coordinates we have chosen. 

All these results can be at once applied to spherical and cylindrical coordi- 
nates, but as this has nothing to do with our present purpose — the exemplification 
of Quaternion methods — we leave the matter here. 

Let us as an example of particular coordinates to which this section forms a 
suitable introduction consider St Venant's Torsion Problem by means of cylindri- 
cal coordinates. 

Saint-Venant's Torsion Problem. 

27. In this problem we consider the equilibrium of a cylinder with any given 
cross-section, subjected to end-couples, but to no bodily forces and no stress on 
the curved surface. 

We shall take r, (p, z as our coordinates, the axis of z being parallel to the 
generating lines of the cylinder. Let A,fi, v be the unit vectors in the directions of 
dr, d(f>, dz respectively and let 

rj = uA + vfi + wv 

as before. 

We shall follow Thomson and Tait's lines of proof — i.e. we shall first find the 
effect of a simple torsion and then add another displacement and so try to get rid 
of stress on the curved surface. 



"[Note added, 1892. In the Phil. Mag. June, 1892, p. 488, there is a mistake in the equation 
just preceding equation 31 and there are two mistakes in equation 31. In the first of these IClD^u- 
n "U/(V — (VJjjw) should be 2I(D^u- n m^v - ^m n w). In equation (31) all the 2's should be dropped.] 
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Holding the section z = fixed let us give the cylinder a small torsion of 
magnitude t, i.e. let us put 

r\ = rzr/u, (40) 

for all points for which rz is small. 

The practical manipulation of such expressions as this is almost always facili- 
tated by considering the general value of Q(V\, 71) where Q is any function linear 
in each of its constituents. Thus in the present case 

G(Vi, ?7i) = t{zQ(A,/u) - zQ(m, X) + rQ(y,n)}. 

[If Q is symmetrical in its constituents, e.g. in the case of stress below this reduces 
to the simple form Q{V\, 771 ) = TrQ(v,/u).] From this we at once see that 77 satisfies 
the equation of internal equilibrium 



for putting Q(a,/3) 



so that both S V77 and V 2 r/ = 0. 

Again the value for Q at once gives us the stress for 

vjcx) = -nScoV . rj - nViStorji - (m - n)a>SVri, 
or urco = -nTr(fjSojv + vS co/j), (41) 

which is a shearing stress nrr on the interfaces perpendicular to fi and v. 

Putting co = the unit normal of the curved surface we have for the surface 
traction 

mco = -nrrv S cofi. 

In the figure let the plane of the paper be z = 0, O the origin, P a point 
on the curved surface and OM the perpendicular from O on the tangent at P. 
Thus -rSco/u = OP cos OPM = PM, PM being reckoned positive or negative 
according as it is in the positive or negative direction of rotation round Oz- Thus 
we see that the surface traction is parallel to Oz and = nrPM. 

M ^\ u 

A 





nV 2 rj + mVSVri = 0, 


a/3 




V?7 = 


= t(2zv - rX) = tV(z 2 - r 2 /2) 
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Hence in the case of a circular cylinder a torsion round the axis satisfies all 
the conditions of our original problem, but this is true in no other case. 

The surface traction at any point on the plane ends necessary to produce this 
strain is mv = riTr/u by equation (41) so that its moment round the origin is 
nrjjr 2 dA, where dA is an element of area and the integral extends over the whole 
cross-section. 

28. Let us now assume a further displacement 

rj = wv, (42) 

where wis a function of r, <p only, and let us try to determine w so that there is 
still internal equilibrium and so that the stress on the curved surface due to w shall 
neutralise the surface traction already considered. 
In the present case 

G(Vi,i7i) = fi(Vw,v). 

Thus SVrj = (since v is perpendicular to Vw) and therefore the equation of 
internal equilibrium gives 

V 2 w = 0. 

Again etw = -nSd>V . r] - nViScorji, 

or nra> = -niyScoVw + VwScov), (43) 

a shear = nTVw on the interfaces perpendicular to Vw and v. Thus putting co for 
the unit normal to the curved surface, the present surface traction will neutralize 
the former if 

StoVw = -TrSco/x = -TScov(Ar) 
= \TSvuN{r 2 ), 
i.e. dw/dn = d(rr 2 /2)/ds, 

where d/dn represents differentiation along the normal outwards and d/ds differ- 
entiation along the positive direction of the bounding curve. 

We leave the problem here to the theory of complex variables and Fourier's 
Theorem. 

Observe however that the surface traction at any point on the plane end is 
equal to mv = nVw by equation (43), and therefore that the total couple is equal 
to -nffrS/j-VwdA = nff(dw/d([>)dA. This leads to the usual expression for 
torsional rigidity. 
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Wires. 

29. In the following general treatment of Wires some of the processes are 
merely Thomson and Tait's translated into their shorter Quaternion forms; others 
are quite different. The two will be easily distinguished by such as are acquainted 
with Thomson and Tait's Nat. Phil. 

The one thing to be specially careful about is the notation and its exact mean- 
ing. This meaning we give at the outset. 

The wires we consider are not necessarily naturally straight but we assume 
some definite straight condition of the wire as the "geometrically normal" condi- 
tion. 

The variable in terms of which we wish to express everything is s the distance 
along the wire from some definite point on it. 

Any element of the wire, since it is only slightly strained, may be assumed to 
have turned as a rigid body from its geometrically normal position. This rotation 
is expressed as usual (Tait's Quaternions, § 354) by the quaternion q; the axis 
of q being the axis of rotation, and the angle of q, half the angle turned through. 

co is taken so that the rate of this turning per unit length of the wire is qcoq' 1 so 
that co is the rate of turning per unit of length when the whole wire is moved as a 
rigid body so as to bring the element under consideration back to its geometrically 
normal position. Of course to is a function of q and its derivative with reference 
to s. This function we shall investigate later. The resolved part of qtoq' 1 parallel 
to the wire is the vector twist and the resolved part perpendicular to the wire is in 
the direction of the binormal and equal to the curvature. In fact co is the vector 
whose coordinates are the k, A, t of Thomson and Tait's Nat. Phil. § 593. When 
we are given q or co for every point we know the strain of the wire completely. 
u>o is defined as the naturally normal value of co, i.e. the value of co when the wire 
is unstressed. 

As usual we take p as the coordinate vector of any point of the wire, p like the 
rest of the functions being considered as a function of s. We shall denote (after 
Tait, Quaternions, Chap, ix.) differentiations with regard to s by dashes. 

We now come to the dynamical symbols. F and M are the force and couple 
respectively exerted across any normal section of the wire on the part of the wire 
which is on the negative side of the section by the part on the positive side. 

Finally let X, L be the external force and couple per unit length exerted upon 
the wire. 
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30. When the wire is strained in any way let us impose a small additional 
strain represented by an increment 6to in to and an increment 8e in the elonga- 
tion at any point. Then the work done on the element ds by the stress-force 
= -deSFp' ds and that done by the stress-couple = -Sqdtoq~ 1 Mds. If (as we 
assume, though the assumption is not justified in some useful applications of the 
general theory of wires) F and M to be of the same order of magnitude the for- 
mer of these expressions can be neglected in comparison with the latter for 6e is 
a quantity small compared with 6to. Now the work done on the element by the 
stress = the increment of the element's potential energy = 5wds where w is some 
function of the strain. Hence 

5w = -S dtoq~ 1 Mq. 
Thus w is a function of to only and 



whence we see that 



6w = -S dtoJJw, 



q l Mq = JJw. (46) 



Notice that M is Thomson and Tait's £, q, £ and q~ l Mq their KLM {Nat. Phil. 
§§594,614). 

31. Now since the strain is small, q^Mq is linear in terms of the strain and 
therefore in terms of to. Hence we see that w is quadratic in terms of to. Let us 
then put 

W = W 2 (tO - OJq, tO - OJq) + W\(tO - to Q ) 

+ a function of temperature only, 

w 2 and w\ being linear and homogeneous in each of their constituents. This is the 
most general quadratic function of to. Now 

J/W = £W 2 (£, 0) - OJq) + £w 2 (oj - to , + £Wi£. 

Putting then to = to and JJw = we get E,W\E, = 0. Operating on the last by 
Scr( ) where a is any vector we see that w\ = 0. Thus putting 

JJ W = (f>(oj - OJq), 
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where cp has the value given by the last equation and is therefore self-conjugate 
we get the two following equations 

w = -S(a> - cl> q )<P(cl) - coq)/2 + wq, (47) 

[as can be seen by a comparison of the last three equations wq being the function 
of the temperature] and 

q~ l Mq = <p(aj - co ). (48) 

When the natural shape of the wire is straight these become 

w = -Sa>(pa>/2 + wq, (49) 

q- l Mq = (pa>, (50) 

and when further the wire is truly uniform <p and Wq are constant along the wire. 

32. Assuming the truth of these restrictions let us conceive a rigid body 
moving about a fixed point which, when placed in a certain position which we 
shall call the normal position, has, if then rotating with any vector angular veloc- 
ity to, a moment of momentum = (pco where <p has the meaning just given. If the 
rigid body be made to take a finite rotation q{ )q~ x and then to move with angular 
velocity qtoq~ l its moment of momentum will be qcpajq' 1 . Now let a point move 
along the wire with unit velocity and let the rigid body so move in unison with 
it that when the moving point reaches the point s the rigid body shall have made 
the rotation represented by q (§ 29 above). Thus by the definition of a> and q its 
angular velocity at any instant is qcoq' 1 and its moment of momentum therefore 
q<p<x>q~ x or M. 

Now consider the equilibrium of the wire when no external force or couple 
acts except at its ends. In this case F is constant throughout and it is easy to see 
(what indeed is a particular case of equation (53) below) that 

M' + Vp'F = 0. (51) 

Interpreting this equation for our rigid body we get as the law which governs its 
motion 

<i(vect. mom. ofmom.)/^ = -Vp'F. 

Thus the rigid body will move as if acted upon by a constant force F at the end 
of the unit vector p' or — since this vector is fixed in the body — as if acted upon 
by a constant force acting through a point fixed in the body. From this kinetic 
analogue of Kirchhoff 's the mathematical problem of the shape of such a wire as 
we are now considering, under the given circumstances, is shewn to be identical 
with the general problem of the pendulum of which the top is a variety. 
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33. We will now give the general equations for any wire under any external 
actions. The comparison of the Quaternion treatment of this with the Cartesian 
as given in Thomson and Tait's Nat. Phil. § 614 seems to me to be all in favour 
of the former. 

The equations of equilibrium of an element ds are with the notation explained 
in § 29 above 

d¥ + Xds = 0, 
VdpF + dM + Lds = 0, 

or dividing by ds 

F' + X = 0, (52) 

Vp'F + M' + L = 0. (53) 

Operating on the last equation by Vp'( ) noting thatp' 2 = -1 and putting Sp'F = 
-T we get 

F=p'T + Vp'{M + L), (54) 

whence by equations (52) and (53) respectively 

X + d{p'T + Vp'(M' + h)}/ds = 0, (55) 

Sp'(M' + L) = 0. (56) 

Now by equation (48) above 

q- l Mq = (p(co - oj ). (48) 

Also by the definition of q 

p = qAq~\ (57) 

where A is some given constant unit vector. Finally as we are about to prove 

to = 2Vq~ l q'. (58) 

It is usual in the Cartesian treatment to leave the problem in the form of equa- 
tions equivalent to the above (55) to (58), 13 scalar equations for the 13 unknown 
scalars of p, T, M, q and to. We can however as we shall directly in the Quater- 
nion treatment quite easily reduce the general problem to one vector and one 
scalar equation involving the four unknown scalars of T and q in terms of which 
all the other unknowns are explicitly given. 
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To prove equation (58)* observe that 

(q + dq)cr(q + dq)~ l = q(cr + Vcocrds)q~ 1 , 

where <x is any vector. The truth of this is seen by noticing that (q+dq)( )(q+dq)~ l 
is the operator that rotates any vector of the element s + ds from its geometrically 
normal position to its strained position. But we can also get to this final position 
by first in the geometrically normal wire making the small strain co ds at the given 
element and then performing the strain of the wire up to the point s. The first 
process is represented on the left of the last equation and the second on the right. 
Thus we get 

q~ l dqcr + crd{q~ l ) . q = Vuxrds, 
or v d(q~ l ) = -q~ l dqq~ l , 

q~ x qcr - crq~ l q = Vcocr, 
whence co = 2V q~ x q' . 

Returning to equations (55) to (58) observe that equations (57) and (58) give 
p' and co as explicit functions of q. Hence by equation (48) 

M = q(j)(2Vq- l q' - co )q~\ (59) 

which gives M also explicitly. Substituting for M and p' in equations (55) and 
(56) we have 

X + — lqAq- l T + VqAq- l (L + —[q <p{2V q~ l q' - co ) q~ l ^\ = 0, (60) 

S qAq- l (L + — [q (f>(2V q~' q' - co )q- 1 ]) = 0, (61) 

which are sufficient equations to determine T and q, whereupon M is given by 
equation (59), co by equation (58) and p' and therefore also p by equation (57). 

Spiral springs can be treated very simply by means of the above equations, 
but we have already devoted sufficient space to this subject. 



This could be deduced from Tait's Quaternions, § 356, equation (2). His e is our qcoq ' and 
his dots our dashes. 



Section IV. 

Electricity and Magnetism. 

Electrostatics 
General Problem. 

34. Merely observing that all the theorems in integration given in the Pre- 
liminary and ivth chapters of Maxwell's treatise on Electricity and Magnetism, 
Part i., are easy particular cases of equations (8) and (9) § 6 above, we pass on to 
the one application of Quaternions that we propose to make in Electrostatics. 

This is to find the most general mechanical results arising from Maxwell's 
theory of Electrostatics, and to see if they can be explained by stress in the di- 
electric. This problem as far as I am aware has not hitherto in all its generality 
been attacked though the most important practical cases have been, as we shall 
see, considered by Maxwell, Helmholtz, Korteweg, Lorberg and Kirchhoff. 

It is necessary first of all to indicate as clearly as possible what I take to be 
Maxwell's theory of Electricity. 

He assumes* all space to be uniformly filled with a certain substance called 
Electricity. Whatever electrical actions take place depend on the continued or past 



*Prof. J. J. Thomson in his paper on Electrical Theories, B.A. Reports, 1885, p. 125, does not 
credit Maxwell with such a definite and circumscribed theory as that described in the text, and 
he is thereby led to find fault with Maxwell's term "Displacement" and points out that there is an 
assumption made with reference to the connection between the true current and this polarisation 
(displacement). He says moreover, "It is rather difficult to see what is meant in Maxwell's Theory 
by the phrase 'Quantity of Electricity.' " None of these remarks are called for if the view I take 
of Maxwell's theory be correct, and these grounds alone I consider sufficient for taking that view. 
The paper of Thomson's here mentioned I shall frequently have to refer to. [Note added, 1892. In 
the text I have given much too rigid a form to Maxwell's theory. What I have called his theory I 
ought rather to have called his analogy. Still I think the present foot-note is in the main just. In my 
opinion it is no more and no less difficult to see what is meant in Maxwell's Theory by "Quantity 
of Electricity" than by "displacement" since the two are connected by perfectly definite equations. 
Of course it is wrong to define "displacement" as "displacement of quantity of electricity" and 
then to define "quantity of electricity" in terms of "displacement" but Maxwell does not seem to 
me even tacitly to do this. Rather he says — the dielectric is polarised; this polarisation can be 
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motion of this substance as an incompressible fluid. If electricity is brought from 
a distance by any means and placed in a given space there must be a displacement 
of the original electricity outwards from that space and the quantity of foreign 
electricity is conveniently measured by the surface integral of that displacement. 

Dielectrics are substances in which this displacement tends to undo itself, so 
to speak, i.e. the original electricity tends to go back to its primitive position. In 
conductors, on the other hand, there is no property distinguishing any imported 
electricity from the original electricity. 

The rate of variation of displacement, whether in dielectrics or conductors, of 
course constitutes an electric current as it is conveniently called. 

We have next to consider a vector at each point of space called the electro- 
motive force, which depends in some way at present undefined on the distribu- 
tion of the displacement in the dielectrics, the distribution of currents whether 
in dielectrics or conductors, and on extra-electrical or semi-electrical action, e.g. 
chemical or mechanical. 

If at any point the electro-motive force be multiplied by a scalar the medium 
at the point remaining (except electrically) unchanged, the current in the case of 
conductors and the displacement in the case of dielectrics is altered in the same 
ratio. In other words the current or the displacement, as the case may be, is 
a linear vector function of the electro-motive force, and the coordinates of the 
linear vector function* at any point depend solely on the state of the medium 
(whether fluid, solid, &c, or again strained or not) at that point. 

To complete the theory we have to explain how the part of the electro-motive 
force which is a function of the distribution of displacement and current depends 
on this distribution. This explanation is obtained by making the assumption that 
the electro-motive force bears to electricity defined as above exactly the same 



represented by a vector D; electrical quantity can be expressed in terms of D; the mathematical 
connections between electrical quantity and D are the same as those between quantity of matter in 
a space and the displacement out of that space made by other matter to make room for the given 
matter; we will impress this useful analogy firmly on our minds by calling D the displacement. 
But I have expressed my present views on the meaning of Maxwell's theory much more fully in 
Phil. Trans. 1892, p. 685.] 

This frequently recurring cumbrous mode of description must be tolerated unless a single 
word can be invented for "a linear vector function of a vector." Might I suggest the term "Hamil- 
tonian?" Thus we should say that the displacement is a Hamiltonian of the electro-motive force, 
the Hamiltonian at any point being a function of the state of the medium. 
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energy relation as ordinary force does to matter, i.e. — 

(wk. dn. on electricity moved during any displacement) 
= (total displacement of elect.) 

x (resolved part of e.m.f. in the direction of displacement). 

In the ivth part of Maxwell's treatise he gives complete investigations of the 
mechanical results flowing from this theory so far as it refers to currents, but he 
has not given the general results in the case of Electrostatics. Nor has he shewn 
satisfactorily, it seems to me, that the ordinary laws of Electrostatics flow from 
his theory. It is these investigations we now propose to make. 

35. Our notation will be as far as possible the same as Maxwell's. Thus for 
the displacement at any point we use D, and for the e.m.f. E. From the connection 
explained in last section between D and E we have 

D = KE/4n, (1) 

where K at any point is some linear vector function depending on the state of the 
medium at the point. If the medium change in any manner not electrical, e.g. by 
means of ordinary strain K will in general also suffer change. 

Let w be the potential energy per unit volume due to the electrical configu- 
ration. Thus if a small increment 5D be given to D at all points, the increment 
JJJ Sw ds in fffw ds, the potential energy of the electrical configuration in any 
space, = work done on the electricity in producing the change, 

fffSwds = -fffSESDds, 

by the relation stated in § 34 existing between E and D. Thus limiting the space 
to the element ds> 

5w = -SESD. (2) 

Now suppose D = nD' so that by equation (1) E = nE' where E', D' are corre- 
sponding e.m.f. and displacement respectively. Thus 

8w = -ndnSE'D' . 

Integrating from n = to n = 1, and finally changing D', E' into D, E we get 

w = -5DE/2. (3) 
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From this we get 

6w = -SE6D/2-SD6E/2, 

so that by equation (2) S E 6D = S D 6E or by equation (1) 

SEK6E = S 6EKE. 

Hence (because E and 6E are quite arbitrary) K is self -conjugate and therefore 
involves only six instead of nine coordinates* . 

In electrostatics the line integral of E round any closed curve must be zero, 
for otherwise making a small conductor coincide with the curve we shall be able 
to maintain a current by § 34, and so (by the same section) constantly do work on 
it (i.e. as a matter of fact create heat) without altering the statical configuration. 
Hence E must have a potential, say v. Thus 

E = -Vv. (4) 

Since in an electrostatic field there is no current in a conductor, E = throughout 
any such conductor and therefore v = const. 

36. The charge in any portion of space is defined as the amount of foreign 
electricity within that space. Thus the charge in any space is the surface integral 
of the displacement outwards. Thus if there be a charge on the element dL of a 
surface in the dielectric this charge = S <TL a D a + S dL h D h where a, b denote the 
two faces of the element (so that <TL U = -<TLb) and in accordance with § 1 above 
dL a points away from the region in which the displacement is D a . Thus <x being 
the surface density 

a- = [SDU JEW, (5) 

where the notation [ ] a+h is used for [ ] a + [ ] h . Similarly if there be finite volume 
density of foreign electricity, i.e. finite volume density of charge in any space, the 
charge = -ffS DdL = -JJJS VD ds, so that if D be the volume density 

D = -SVD. (6) 

[The reason for having +S <TL a D a before and -S dL D here is that in the former 
case we were considering a charge outside the region where D a is considered — 
between the regions a and b in fact — whereas in the latter case we are considering 



*We see from this that D = eVw or E = dVw according as w is looked upon as a function of 
EorD. 
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the charge inside the region where D is considered. The same explanation applies 
to the sign of -SDU dL for the surface of a conductor given below.] 

In conductors, as we saw in § 34, the displacement has virtually no meaning 
(except when it is changing and so the phenomenon of a current takes place) 
for the foreign electricity and the original electricity are not to be distinguished. 
Not so however with the surface of the dielectric in contact with the conductor. 
We may therefore regard the electricity within the body of the conductor as the 
original electricity so that the charge is entirely at the surface. Thus the surface 
density will be -SDU dL where dl, points away from the conductor and D is the 
displacement in the dielectric. This may be regarded as a particular case of equa- 
tion (5) D being in accordance with what we have just said considered as zero in 
the conductor. 

37. All the volume integrals with which we now have to deal may be con- 
sidered either to refer to the whole of space or only to the dielectrics, as the con- 
ductors (except at their surfaces) in all cases contribute nothing. The boundary of 
space will be considered as a surface at infinity and all surfaces where either D 
or E is discontinuous. 

Putting W = fffw ds we have already found one expression for W, viz. 

2W = -fffSDEd$. 

We now give another. By equation (4) 

2W = fffSDVvds 

= ffvSdUD- fffvSVDds, 

by equation (9) § 6 above. Thus by equations (5) and (6) § 36, 

2 W = ffvcr ds + fffvD d%, (7) 

where ds is put, as it frequently will be, for an element of surface, i.e. T dL. The 
value of W which we shall use* is obtained by combining these two, viz. 

W = ffvcrds + fffvDds+^fffSDEds. (8) 



This is for the general case following the example of Helmholtz in the particular case when K 
reduces to a single scalar. See Wiss. Abh. vol. i. equation (2d), p. 805. The method adopted in the 
following investigation is also similar to his. 
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So far we have merely been shewing that all the above results of Maxwell's 
flow from what in § 34 has been described as his theory. We now proceed to the 
actual problem in hand which is proved from these results however they may be 
obtained. I may remark that some such investigation as the above seems to me 
necessary to make the logic of Maxwell's treatise complete. 

38. Suppose now that W is the potential energy of some dynamical system 
extending throughout space. Let us give to every point of space a small displace- 
ment 6n vanishing at infinity and find the consequent increment 6W in W. If this 
can be put in the form 

5W = -fffS 6t} X 0Vj ds, (9) 

we shall have the following expression for F the force per unit volume due to the 
system 

F = <f>A, (10) 

and the following expression for ¥ s the force per unit surface at any surface of 
discontinuity in (p 

¥ s = -[cf>UdY] a+b , (11) 

the notation being the same as in equation (5) § 36. 

Moreover if $ be self-conjugate the forces both throughout the volume and 
at surfaces of discontinuity are producible by the stress 4> as can be seen by § 1 3 
above. [Compare all these statements with § 14 above.] 

For proof, we have by equation (9) § 6 

5W = -fffS 6^ 0Vi ds = -ffS dr/(/)di: + fffS 6rj fa Vi d% 

where of course the element d~L is taken twice, i.e. once for each face. But 

SW = -(work done by the system F, F s of forces) 
= ffS 6rjF s ds + fffS 6t]F ds, 

where the element ds is taken only once. Equating the coefficients of the arbitrary 
vector 5r] for each point of space we get the required equations (10) and (11). 

39. We must then put SW where W is given by equation (8) in the form 
given in equation (9). 

We must first define 8 when applied to a function of the position of a point. 
Suppose by means of the small displacement 6n any point P moves to P' . Then Q 
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being the value at P, before the displacement 6rj, of a function of the position of 
a point, Q + 6Q is defined as the value of the function at P' after the displace- 
ment. Thus even in the neighbourhood of a surface of discontinuity 6Q is a small 
quantity of the same order as 5i]. 

Now the charge within any space, that is the quantity of foreign electricity 
within that space will not be altered by the strain. 

6(Dds) = 0, 6((rds) = 0. (12) 

To find SV we have 

S(dp + 6dp)(V + <5V). = S dpV. 

or, since 6dp = -S dpV . 8rj, 

S dpdV = S dpViS cfyiV, 

whence 5V = V : S 5^ V. (13) 

The part of 6W depending on the first two terms of equation (8) is by equa- 
tion (12) 

Jffdv Dds + ffdv a ds 

= -fffSv SVDds + ffSv S dZD [by equations (5) and (6) § 36] 

= fffS D V 5v ds [by equation (9) § 6] . 

Noticing that 6ds = -dsSVdr] and that 4ndD = K6E + 6KE we see that the last 
term in equation (8) contributes 

fffS D SE ds - ijjjf S DES V 6t] ds + (&ny l fffS E 6K E ds. 

Combining the last result with the first term of this we get 

fffS D( V 6v - <5(Vv)) ds = - fffS D 6V . v ds 
= -JffSDV 1 S6r] l Vvds 
= fffSDViSESriids. 

Thus we have 

6W = -fffS dmijViSDE -ESV^ds + (8nT l fffSESKEds. (14) 
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40. Now the increment 6K in K is caused by two things viz. the mere rota- 
tion of the body and the change of shape of the body. Let us call these parts 6K r 
and 5K S respectively. 

First consider 8K r . Suppose the rotation is e so that any vector which was to 
becomes thereby to + Veto. Thus the result of operating on to + Veto by K + 5K r 
is the same as first operating on to by K and then rotating. In symbols 

(K + 6K r ){co + Veto) = Kto + VeKto, 

whence 8K r co = VeKco - K(Veco). 

Thus SE6K r E = SEeKE - SEK(VeE) 

= 2SeKEE = SnS eDE, 

whence giving 6 its value WSrj/2, 

(S7T)- l SESK r E = SV 6ri VDE/2. 

Substituting in equation (14) 

6W = -JffS 577i{yV x DE/2 -ESVjD} ds + (8^) _1 fffSE6K s Eds, 

or 

6W = ±fffS 5rn VD V X E ds + (Sny 1 fffSE6K s E ds. (15) 

It only remains to consider 6K S . K is a function of the pure strain of the 
medium and 5K S is the increment in K due to the increment in pure strain owing 
to 6rj. Calling this increment of pure strain 5ij/ so that by equation (3) § 12 above 

Sij/co = -Srji ScoVi/2 - ViScodr/i/2, 

we have 

5K S = -S dijf ^a 2 ^ . K 2 , (16) 

by equation (7) § 5 above. This gives by equation (5) § 3 

6K, = -S6t 1 ^Q 2 V 1 .K 2 , (17) 

SE6K S E = -S 6r/ llp a2^iSEK 2 E. (18) 

Now by equations (1) and (3) § 35, 

8nw = -SEKE, 
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so that w is a function of the independent variables E, if/ (because K is a function 
of if/). Therefore 

(Sny l SE6K s E = S6tj h qwV 1 . (19) 

This equation might have been deduced at once thus 

(8n)- 1 SE6K,E = S Sif/^Qw$ = S Srj^QwVu 

but equation (17) is itself of importance so the above proof is preferable. 
Thus finally from equation (15) we get 

8W = J3j5(J^i(VDViE/2 + ^awVi)ds. (20) 

We therefore have for in equations (10) and (11) § 38, 

<t>u) = -VDa)E/2 -^awu). (21) 

This is a self-conjugate function so that as we saw in § 38 it is a stress which 
serves to explain forces both throughout the volume of the dielectric and over any 
surfaces of discontinuity in D or E*. 

The force in particular cases. 

41. Let us first consider that part -VDAE/2 of the force (equations (10) 
§ 38 and (21) § 40) which does not depend on the variation of K with the shape 
of the body. 

Suppose our dielectric is homogeneous and electrically isotropic so that K is 
a simple constant scalar. In this case 

4nD = -KVv, (22) 

by equations (1) and (4) § 35. Therefore by equations (5) and (6) § 36, 

4ttD = KV 2 v, (23) 

4no- = -K[S U dY. Vv] a+b . (24) 

From these we at once get by the theory of potential that 

Kv = fffuD ds + ffuo- ds. (25) 



*As far as I am aware nobody has hitherto attempted to find the electrical forces much less 
the stress except in the case when D is parallel to E i.e. the dielectric is electrically isotropic 
when unstrained. The particular results contained in § 45 below have been obtained by Korteweg, 
Lorberg and Kirchhoff as is stated in Prof. J. J. Thomson's paper (p. 155) referred to in § 34. 
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From this we know by the theory of potential that at the surface where the charge <x 
resides Vv is discontinuous only with regard to its normal component and at all 
other points is continuous. Thus 

Vv a = Vv b + xUdL a 

and by equation (24) x = Ana/K so that 

(4nT 1 K(Vv a -Vv h ) = crUdi: a , 
whence (4n)~ l KU dL a (Vv a - Vv fo ) = -cr. (26) 

Now the force F per unit volume is 

-VDAE/2 = -KVVvAVv/8n 
= -KV 2 vVv/4n, 
or F = -DVv, (27) 

and the force per unit surface F s . is by equation (11) § 38, 

\[VDUdYE] a+b = (8n)- l K[VVvUdi:Vv] a+b 

= (SnT'KViVvaUdZa Vv a - Vv h U dL a Vv,) 
= {%ny l KV{Vv a + Vv h )U dL a (Vv a - Vv,), 

whence by equation (26) 

F, = -±(r[VvW (28) 

Thus we see that Maxwell's theory as given in § 34 above reduces to the 
ordinary theory when K is a single scalar. In fact two particles containing charges 
e e' apparently repel one another with a force ee' /Kr 2 where r is the distance 
between them, for by equations (25), (27) and (28) the force in any charged body 
is that due to a field of potential v given by 

Kv = Y.ue. (29) 

42. If the medium when strained remain electrically isotropic ^(jK as well 
as K must be a simple scalar. Thus with Thomson and Tait's notation for strain, 
which makes the coordinates of if/, e, f, g, a/2, b/2, c/2 we have 

^Q^ = dK/de = dK/df = dK/dg 
dK/da = dK/db = dK/dc = 0. 
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Therefore K is a function of e + / + g only, i.e. of the density (m) of the medium. 
Thus because 

de + df + dg = -dm/m = -d log m 

we get ^qK = -dK/d log m = -k suppose. Hence 

^aw = (SnT^EidK/d log m)E = kE 2 /8n. 

Thus the force -^QwA [equations (21) §40 and (10) §38] resulting from the 
change of K with pure strain is in the case we are now considering 

- kVE 2 /8n, (30) 

and is* therefore, according as k is positive or negative, in the direction of or that 
opposite to that of the most rapid increase of the square of the electro-motive 
force. Thus even in the case of a fluid dielectric which has no internal charge 
but which forms part of a non-uniform field of (electro-motive) force the surfaces 
of equal pressure and therefore the free surface will if originally plane no longer 
remain so. 

Nature of the Stress. 

43. We have seen that the stress which serves to explain the electrostatic 
forces is that given by equation (21) § 40, viz. 

cpco = -VDcoE/2- l! ,awco. (21) 

Let us first consider the part -VDcoE/2 which does not depend on the varia- 
tion of K. Putting co first = £/D and then = C/E we get 

</>UD = TDTE . f/E/2, 
ciUE = TDTE . t/D/2. 

Therefore putting co first = any multiple of f/D + UE and then = any multiple of 
f/D - UE we get 

cpco = TDTEco/2, 
(poo = -TDTEco/2. 



This is the same result as Helmholtz's on the same assumption Wiss. Abh. i. p. 798. 
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Lastly, since 

-VDcoE = coSDE - DScoE - EScoD, 

we see that if we put to = any multiple of VDE 

<fito = toSDE/2 = -cow. 

Thus we see that the stress now considered is a tension along one of the bisec- 
tors of D and E (the bisector of the positive directions or the negative directions 
of both) = TDTE/2, an equal pressure along the other bisector and a pressure 
= w perpendicular to both these directions. When D is parallel to E this at once 
reduces to Maxwell's case, viz. a tension in the direction of E and a pressure in 
all directions at right angles each = w. 

44. We have now to consider the other part of the stress, viz. 

<pto = -^(jwto, 
or cpaj = ^QiOjSEKiE/in. (31) 

If we assume that K is a function of the density (m) of the medium only we shall 
have 

dK/de = dK/df = dK/dg = -dK/dlogm = -k, 

say, and 

dK/da = dK/db = dK/dc = 0, 

as in § 42. Here however k is not in general a mere scalar but a self-conjugate 
linear vector function. We have then in this case 

(pto = -coSEkE/Sn, 

which is a hydrostatic pressure or an equal tension in all directions according as 
SEkE is positive or negative. In this case the 36 coordinates of ^Qj to . K\ reduce 
to the 6 of k for each point of space. 

A more general assumption is that 6K S (§ 40) depends only on the elongations 
in the directions of the principal axes of K. Taking i, j, k as unit vectors in these 
directions we again have 

dK/da = dK/db = dK/dc = 0, 
and thus (pi = iSE(dK/de)E/8n, 

and similarly for j and k, so that the principal axes of the stress now considered 
are the principal axes of K. 



§ 46. ] ELECTRICITY AND MAGNETISM. 67 

45*. The most natural simple assumption for solid dielectrics seems to me to 
be that the medium is electrically isotropic before strain, and also isotropic with 
regard to the strain in the sense that if the strain be, so to speak, merely rotated, 
6K S will suffer exactly the same rotation. We may treat this problem exactly as 
we did (§ 20) that of stress in terms of strain for an isotropic solid. Thus splitting 
up 6if/ into its principal elongations, i.e. putting 

dif/co = -UdeiSiu), 

we shall get, as in § 20, 

6K S co = -(a - /?)£ 6ei S ico + j3ufL 5e 
= (a - 0)6ifr co - ficoS £ 6ifr {. 

But 5K s co = -S 5if/ ^di£ . Kico by equation (16) § 40, so that from equation (31) 

-BnS 6iff £<p£ = (a-/3)SE5i//E- /3E 2 S £ Sift (, 

whence we see by § 4 above that 

cpco = {(a-p)ESEco +{3E 2 co}/8n, (32) 

which consists of a pressure in the directions of the lines of force = -aE 2 /8n and 
another pressure in all directions at right angles = -f3E 2 /8n. 

Magnetism 

Magnetic potential, force, induction. 

46. We now go on to the ordinary theory of magnetism; and here we shall 
merely follow Maxwell in his General Theory, so as to give an opportunity of 
comparing Quaternion proofs with Cartesian, as we have already done in Elastic- 
ity. 

We shall not consider in detail the effect of one small magnet upon another, 
as this has already been done by Tait. In connection with this I am content to 
remark that I think the treatment of this problem can be made somewhat simpler 
than Tait's by means of potential. 

Suppose we have a pole —m at O and a pole +m at O' where 00' is small. 
Calling the vector from O to O' 00' , let us call the vector mOO' p, so that p is 



*For references to former proofs of this see foot-note to § 40 above. 
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the vector magnetic moment of the magnet. The potential of -m at any point P 
is -mu, where u as usual = PO l . Similarly the potential of +m is mu' , where 
u' = PO ' 1 . Therefore the potential of the magnet 



= m(u' - u) = mS OO'Vu, 

where of course P is the variable point implied by V. Thus the potential of a small 
magnet /i at any point = S/uVu. 

Hence the potential of any magnet whose magnetic moment per unit volume 
at any point is I is 

Q = SVfffulds = -fffSlVuds, (33) 

according to the convention of § 9 above. By equation (9) § 6 this may be put 

D. = -JfuS I dH + fffuS VI ds, (34) 

which shews that we may consider it due to a volume density S VI and a surface 
density -5 If/ dl,* of magnetic matter, the surface density occurring wherever 
there is discontinuity in I. 

By again considering the poles m and -m of the small magnet // we see that 
its potential energy when placed in a field of magnetic potential Q. is -S/uVQ., 
whence just as we obtained equation (33) we now see that the potential energy 
(W) of any magnet in such a field is 

W = -fffSF7£lds, (35) 

or W = -ffSlSIdL + fffnSVIds, (36) 

by equation (9) § 6 above, so that the potential energy is just the same as it would 
be for the imaginary distribution of magnetic matter. 

47. The force (H) on a unit magnetic pole at any point external to the magnet 
is given by 

H = -VQ = -VSVfffulds = -V 2 fffulds + VVVfffuIds 
= -4ttI + VfffVIVuds = -4nl + VA, 

where 

A = fffVIVu ds. (37) 



"[Note added, 1892. More generally and better -[SIU dL] a+b .] 
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Thus we see that for all external points H = VA, so that A is called the vector 
magnetic potential. [It is to be observed that since VA = H + 4^1 = a vector, 
SVA = 0.] VA is called the magnetic induction and for it we use the single 
symbol B so that 

B = VA. (38) 

Thus SVB = SV 2 A = 0, (39) 

and also by the equation for H just given 

B = H + 4ttI. (40) 

This is not the way in which Maxwell defines the magnetic force and induc- 
tion, but he shews quite simply (Elect, and Mag. §§ 398-9) that his definition and 
the present one are identical. This can be shewn as easily without analysis at all. 

48. Where I is discontinuous both H and B are also discontinuous. From 
the surface density view we gave in equation (34) we see that, just as we have the 
expression, given in § 41 for Vv a - Vv^, so now 

U h -U a = -4nU dZ a [SW dY\ a+b , (41) 

so that the discontinuity in H is entirely normal to the surface of discontinuity. 
Further from this equation we have 

S(H b --a a )dZ a = 4n[SIdS\a + b, 

i.e. [S(U + 4nl)di:] a+h = 0, 

or [SBdX] a+b = (42) 

so that the discontinuity in B is entirely tangential. 

From this equation we see that for any closed surface whatever whether it 
include surfaces of discontinuity in I or not 

ffS d2B = 0. 

For adding these surfaces to the boundary of the inclosed space, in accordance 
with § 7 above, we see by equation (42) that they contribute zero to the surface 
integral; but the total surface integral is by equation (9) § 6 JJJSVBds = by 
equation (39). 
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Magnetic Solenoids and Shells. 

49. A magnet is said to be solenoidal if the imaginary magnetic matter of 
equation (34) is entirely on the surface. Thus for a solenoidal distribution 

SVI = 0. (43) 

In this case the potential is by equation (34) 

Q = -ffuSldZ. (44) 

50. A simple magnetic shell is defined as a sheet magnetised everywhere 
normally to itself and such that, at any point, the magnetic moment per unit sur- 
face is a constant called the strength of the shell. 

Calling the strength cp we have for the potential energy at any point by equa- 
tion (33) 

a = -cf) ffs dL Vw. (45) 

Now -S dL Vu is the solid angle subtended by the element dL at the point con- 
sidered, so that 

D. = (p x (solid angle subtended by shell at point). (46) 

Thus if P be a point on the positive side of the shell and P' a point infinitely 
near P but on the negative side 

Potential at P- Potential at P' = Ancp, 

or what comes to the same thing 



I 



SHdp = 47T<p. 



This integral may be taken along any path, e.g. along a path which nowhere cuts 
the shell. The same integral is true if H be the magnetic force due to a whole field 
of which the shell is only one of several causes, for the part contributed by the 
rest of the field is zero on account of the infinite proximity of P and P' . For future 
use in electro-magnetism observe that this statement cannot be made if for H in 
the integral be substituted B. 
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51. The condition that any magnet can be divided up into such shells is at 
once seen to be that I can be put in the form 

I = V0, (47) 

where (p is some scalar. 

In this case the potential is by equation (33) 

Q. = -fffsiVuds = -JffSV(f>Vuds, 

or by equation (9) § 6 

Q = -ff<f>S dX Vw + fff(f>V 2 u ds, 
i.e. Q. = -ff<f>S dXVu + 4n<p. (48) 

Remarking that the solid angle again occurs here it is needless to interpret the 
equation further. By equation (37) we have for the vector potential 

A = fffVIVuds = fffWcpVuds, 

or by equation (9) § 6 

A = fffiVdHVu. (49) 

52. The potential energy of a magnetic shell of strength cp placed in a field 
of potential Q, is of importance. We see by equation (35) that it is 

W = -<f)ffS cE va 

If then the magnets which cause Q. do not cut the shell anywhere, so that -VQ, = 
VA, we shall have 

W = (pffS dJ. VA = ffS dIB, (50) 

or W = <f> I S dp A, (51) 

by equation (8) § 6. 

Suppose now that A is caused by another shell of strength cp' . Then by equa- 
tion (49) 

A = (P'ffVdZ'Vu = <p' I udp', 

by equation (8) § 6. Thus finally the potential energy M of these two shells is 
given by 

M = <fxp'ffuS dp dp'. (52) 
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53. The general theory of induced magnetism when once the proposition 
(given in equation (42) § 48) that [S dl,B] a+h is zero is established, is much the 
same whether treated by Quaternion or Cartesian notation. We shall therefore not 
enter into this part of the subject. 

Electro-magnetism 

General theory. 

54. We now propose to prove the geometrical theorems connected with 
Maxwell's general theory of Electro-magnetism by means of Quaternions. 

We assume the dynamical results of Chaps. V., VI. and VII., and the first six 
paragraphs of Chap. VIII. of the fourth part of his treatise. 

These assumptions amount to the following. Connected with any closed curve 
in an electro-magnetic field there is a function 



-/ 



SXdp, (53) 



where A is some vector function at every point of the field. The function p has 
the following properties. If any circuit be made to coincide with the curve the 
generalised force acting upon the electricity in the circuit is 

E = -p. (54) 

Again, if there be a current of electricity y flowing round this circuit, the gen- 
eralised force X, corresponding to any coordinate x of the position of the circuit 
due to the field acting upon the conductor, is 

X = ydp/dx. (55) 

55. The first thing to be noticed is that p can be transformed into a surface 
integral by equation (8) § 6 above. 

Thus p = -ffSBdL, (56) 

where B = VVA, (57) 

so that SVB = 0. 

Next we see by the fundamental connection (§ 34 above) between the e.m.f. E 
and electricity, that E must equal the line integral of E round the circuit, or 

E = -JSEdp. (58) 

We are now in a position to find E in terms of A and B, i.e. of A. 
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56. The rate of variation of p is due to two causes, viz. the variation of the 
field (A) and the motion of the circuit (p). In the time 5t then there will be an 
increment 5X in A and an increment 8(TL in <TL to be considered. Thus 

p6t= -JSSAdp-ffSBSdL. 

[This amounts to assuming that J SSAdp = JfSdBdl,, which of course is true 
by equation (8) § 6.] Now when the circuit changes slightly we may suppose the 
surface over which the new integral extends to coincide with the original surface 
and a small strip at the boundary traced out by the motion (p5t) of the boundary. 
Thus ddl, is zero everywhere except at the boundary and there it 

= V(p6t) dp, 

so that p5t = - J S SAdp + 5t J SpBdp, 

whence dividing by 6t 

p = J S (-A + VpB) dp, (59) 

but by equations (54) and (58) p = J SEdp. Thus 

E = -A + VpB - VtA, (60) 

where if/ is a scalar and —Viff is put as the most general vector whose line integral 
round any closed curve is zero. 

57. We now come to the mechanical forces exerted on an element through 
which a current C per unit volume flows. 

We see by equation (55) that the work done by the mechanical forces on any 
circuit through which a current of magnitude y flows in any small displacement 
of the circuit equals yx the increment in p caused by the displacement. Give then 
to each element dp of the circuit an arbitrary small displacement 8p and let F' be 
the mechanical force exerted by the field upon the element. Thus as in last section 

-fSF'dp = y5p = -rffSBSdZ = -y J S 8pdpB. 

Thus the force F' on the element dp is yV dp B. But we may suppose this element 
to be an element ds of volume through which the current C flows. Thus if for y dp 
we write C ds>, and for F', F ds, where F is the force per unit volume exerted by 
the field, we get 

F = VCB. (61) 
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58. So far we have been able to go by considering the electric field as a me- 
chanical system, but to go further (as Maxwell points out) and find how B or A 
depends on the distribution of current and displacement in the field we must ap- 
peal to experiment. It has been shewn by experiment that a small circuit produces 
exactly the same mechanical effects on magnets as would a small magnet, at the 
same point as the circuit, placed with its positive pole pointing in the direction of 
the positive normal to the plane of the circuit when the positive direction round 
the circuit is taken as that of the current*. Moreover the magnetic moment of the 
magnet which must be placed there is proportional to the strength of the current x 
the area of the circuit. Further, the effect of this circuit upon other such small cir- 
cuits is the same as the mutual effects of corresponding magnets. We have now 
only to consider a finite circuit as split up in the usual way into a number of el- 
ementary circuits to see that a finite circuit will act upon magnets or upon other 
circuits exactly like a magnetic shell of strength proportional to the strength of the 
current and boundary coinciding with that of the circuit. The unit current in the 
electro-magnetic system is so taken as to make this proportionality an equality. 

The one difference between the circuit and the magnetic shell is that there is 
no discontinuity in the magnetic potential in going round the circuit, so that by 
§ 50 above the line integral of H round the circuit will be 4nx the strength of the 
current. In symbols 

JSHdp = 4nffSCdi: 

for any curve, so that by equation (8) § 6 

4ttC = VVH, (62) 

whence SVC = 0, (63) 

which of course is a direct result of our original assumption that electricity moves 
like an incompressible fluid. Maxwell tacitly assumes this by making the assump- 
tion that only one coordinate is required to express the motion of electricity in a 
circuit. 

59. We are now in a position to identify the B we are now using with the 
magnetic induction for which we have already used the same symbol. 



*See § 1 above for the convention with respect to the relation between the positive side of a 
surface and the positive direction round its boundary. Hitherto we have had no reason for choosing 
either the right-handed or the left-handed screw as the type of positive and negative rotation. But 
to make the statement in the text correct we must take the former. 
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We see by equation (50) § 52 that the mechanical force on the shell corre- 
sponding to any coordinate x is 

-<f>dffS dZB'/dx, 

where B' is the magnetic induction; and by equations (53) and (55) that the force 
on the corresponding electric circuit is 

-cpdffS dZB/dx, 

therefore B = B' wherever there is no magnetism. And where there is magnetism 
B is not = H for S VB = 0, as we have seen. Thus B = B' at all points. In other 
words the two vectors are identical and we are justified in using the same symbol 
for the two. 

This practically ends the general theory of electro-magnetism. We content 
ourselves with one more application of Quaternions in this subject. We give it 
because it exhibits in a striking manner the advantages of Quaternion methods. 

Electro-magnetic phenomena explained by Stress. 

60. In § 46 equation (35) we have seen that the potential energy of a mag- 
netic element = Ids in a magnetic field is SIHds when H has a potential. 
Maxwell assumes that the same expression is true whether H have a potential 
or not. Assuming this point* with him we can find the force and couple acting on 
the medium and a stress which will produce that force and couple. The force and 
couple due to the magnetism of an element is obtained by giving the element an 
arbitrary translation and rotation and assuming that the work done by this force 
and couple = the decrement in the potential energy of the element. Thus the force 
per unit volume is -ViSIHi for the decrement in the potential energy due to a 
small translation Sp is S Sp ViSlHi. Similarly the couple M is given by 

M = VIH (64) 



*I do not defend the legitimacy of this assumption. It seems to me bold to assume that a mag- 
net possesses any such thing as potential energy in a field which has no potential. If we assume 
H and its derivatives to be continuous throughout our typical element ds of volume containing 
a great number of molecules (both material and magnetic) the force on a magnetic molecule fi 
consisting of two poles is -SfiV . H and the force per unit volume -SIV . H, which is only identi- 
cal with the expression — ViSIHi obtained below when H has a potential. With this expression a 
stress cannot be found that produces the force. If, however, H and its derivatives be not assumed 
continuous in this manner the force on the magnet jj is quite indeterminate whether the magnetic 
pole or the molecular current view of magnetism be taken, unless it be specified in what way the 
poles and currents are distributed in the element of volume. 
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for the decrement -S M 6co in the potential energy due to a small rotation 5a> is 
-S V 5co I . H = —S 6co VIH. The total force F per unit volume is the sum of that 
just given and that given by equation (61), so that 

F = yCB-V^fflj. 

Therefore by equation (62) 

4ttF = VVVH . B - 4ttViSIHi 

= -HiSViB + ViSBHj - 4ttViSIHi 
= -H 1 5V 1 B + V 1 ,SHH 1 . 

Now SVB = so that H 1( S VjB = HSAB, and again 

V^HHi = V(H 2 )/2, 
so that F = <pA, (65) 

where Bncpo) = -2HScoB + coU 2 . (66) 

From this we get 

inVCrt = 2VBH = 8ttV1H, 
so that M = V^rt- ( 67 ) 

From these two results (equations (65) and (67)) we see by § 13 above that 
the stress cp will produce all the mechanical effects of the field. 

This stress, as can be seen by giving co the required values in equation (66), 
is one of pressure -H 2 /87r in all directions at right angles to B and of tension 
-5H(2B - H)/8tt in the direction of H. When there is no magnetism H = B so 
that this pressure and tension become equal and their directions at right angles to 
and along B respectively. In fact we then have 

8n<f>co = -BwB. (68) 



Section V. 
Hydrodynamics. 



61. In the applications I am about to make in this I have practically nothing 
new to shew except the utility of Quaternion methods in the general theory of 
Hydrodynamics in all its parts. 

I therefore take a treatise (Greenhill's article in the Encyc. Brit, on this sub- 
ject) and work out the general theory on the lines of the treatise. This is more 
necessary than at first sight it would seem, for I believe mathematicians who 
have studied Quaternions are under the impression that the method does not lend 
itself conveniently to the establishment and treatment of such equations as the La- 
grangian and those of Cauchy's integrals. With our meaning of V, however, the 
Quaternion treatment of these equations is as much simpler than the Cartesian as 
in the case of the Eulerian equations. 

Down to equation (13) below the subject has already been handled by Prof. 
Hicks* in his Quaternion treatment of Strains and Fluid motion (Quart. Journ. 
Math. xiv. [1877] p. 271). I do not hesitate to go over the ground again as my 
methods are different from his. 

Notation. 

62. For the vector velocity at any point we shall use <x, for the density D, for 
the force per unit mass F, for its potential, when it has one, v, and for the pres- 
sure p. For time-flux which follows a particle we shall use d/dt or the Newtonian 
dot, and for that which refers to a fixed point of space d/dtK 

Thus 

d/dt = d/dt -So-V. (1) 



"[Note added, 1892. Prof. Tait's name ought to be added to Prof. Hicks's.] 
1 [Note added, 1892. I am aware that this is contrary to the usual English custom, but that 
custom — of interchanging the meanings of d/dt and d/dt as given in the text — seems to me out 
of harmony with the meaning attached to d in other branches of Mathematics. At any rate I have 
respectable fellow-sinners, e.g. KirchhofFin his Mechanik, zweite Vorlesung, et seq.] 
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This equation is given on p. 446 of Greenhill's article already mentioned. In 
future we shall refer to this article simply as "Greenhill's article." 

Euler 's Equations. 

63. To find the equation of continuity, with Greenhill, we merely express 
symbolically that the rate of increase of the mass of the fluid in any space equals 
the rate at which it is flowing through the boundary. Thus M being the mass in 
any space, 

dM/dt = ffDScrdZ. (2) 

This is Greenhill's equation (1) p. 445. By it and equation (9) § 6 above we have 

dM/dt = ffJSV(Dcr)ds, 
whence reducing the volume to the element ds, 

dD/dt = SV(Dcr). (3) 

This is Greenhill's equation (2). Thus by equation (1) of last section, 

D/D = SVo-. (4) 

64. To obtain Euler's equations of motion we express that the vector sum of 
the impressed forces for any volume equals the vector sum of the bodily forces 
plus the vector sum of the pressures on the surface. Thus 

SffD&d* = JffDFd*-ffpdL 
= fff(DF-Vp)ds, 

by equation (9) § 6 above. Applying this to the element d$ and dividing by D ds 
we have 

& = F - Vp/D, (5) 

or by equation (1) § 62 

da/dt - S o-V . a = F - Vp/D. (6) 

This is Greenhill's equations (4) (5) (6). 

If F(p, t) = (F a scalar) be the equation of a surface always containing the 
same particles F = 0, or by equation (1) 

dF/dt-So-VF = 0. (7) 

This is Greenhill's equation (7). 
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65. Let us now put 

fdp/D = P. (8) 

This of course assumes that D is a function of p only, which is not always the 
case, for instance in a gas where diffusion of heat is taking place. If F have a 
potential v, F = -Vv. Thus equation (5) becomes 

cr = -V(v + P), (9) 

and equation (6) 

dcr/dt-ScrV . cr = -V(v + P). (10) 

Now 

ScrVr . o-i = Vo-WiO"! + ViSo-o-i = VcrWa- + V(o- 2 )/2. 

Thus equation (10) becomes 

da/dt + IVea + VR = 0, (11) 

where the scalar R is put for P + v - cr 2 /2 and the vector e for Wcr/2. This is 
Greenhill's equations (8) (9) (10). 

If 6 = 0, i.e. W<x = 0, we have cr = V(p, whence our equation becomes 

V(d(/)/dt + R) = 0, 
so that d(f>/dt + R = H, (12) 

where H is a function of t only. In next section we shall obtain in the case of 
an infinite fluid a generalisation of this which I believe has not hitherto been 
obtained. Here we have made the assumption that if Wcr = at one epoch it 
will be so always. This we shall prove later. 

66. Greenhill next considers the case of steady motion. In this case dcr/dt = 
0, so that equation (11) becomes 

2Vecr + VR = 0, (13) 

and therefore the surface R = const, contains both vortices and stream-lines and 
the relation dR/dn = 2qto sin 6 given on p. 446 of Greenhill's article is the natural 
interpretation of our equation. 

So far we have been going over much the same ground as Hicks, but now we 
enter upon applications of Quaternions that I think have not been made before. 
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67. Greenhill next considers rotating axes and finds the form of the equa- 
tions of motion when referred to these. Let <x be the velocity referred to them; so 
that if the axes have at any time made the rotation q( )q~ l the real velocity will 
be qcrq' 1 . Thus*, as always with rotating axes, if a be any vector function of a 
particle the rate of increase of a in space is 

a + Vcoa, 

where co is the angular velocity referred to the same system of rotating axes. Thus 
we see that the acceleration of a particle = & + Vcocr. The velocity cr = p' + Vcop', 
where p' is the vector coordinate referred to our present axes (i.e. in the notation 
of the footnote p' = qpq' 1 ). Thus our equation of motion is now 

& + Vcocr = F - p Np/D, 

but now 

d/dt = d/dt - Sp' p N. = d/dt -S(cr- cop') p N, 

whence changing p' into p we have 

dcr/dt -S(lt- cop)V . cr + Vcocr = F - Vp/D, (14) 

which is the equation Greenhill gives on p. 446. 



This is quoted as a known result because it occurs generally in the subject of Rigid Dynamics. 
No Quaternion proof, however as far as I am aware, has been given. We therefore give one here. 
What is meant by rotating axes may be thus explained. — Instead of choosing as our coordinates 
the vectors a', /?'... which occur in any problem, we take others a, f3 . . . such that a' = qaq^ 1 , 
j3' = q[Sq~ l ■ ■ ■ q( )q~ l may be called the integral rotation of the axes. Thus if we say that 
the vector angular velocity of these axes themselves is a> we mean that the real angular velocity 
is qu>q~ l , so that, as can be seen by putting in equation (a) below a = const., or as in equation (58) 
§ 33 above, a> = 2Vq~ l q. (This maybe established also by Tait's Quaternions, § 356, equation (2), 
from which qcoq~ l = 2Vqq~ l orw = 2Vq~ l q). Again, when we say that the rate of increase of a 
in space is t, we mean that a' = qrq~ l or t = q~ l . d(qaq~ l )/dt . q, or 

t = a + 2VVq~ q . a, (a) 

or t = a + Vcoa. (b) 

This could have been proved with fewer symbols and more explanation, but the above seems to 
me the most characteristic Quaternion proof. We might have started with not quite so general an 
explanation of reference to rotating axes and so refrained from introducing the integral rotation, 
and therefore also q. 
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The Lagrangian Equations* . 

68. We now consider the history of a single particle and for this we require 
different notation. 

We consider the vector coordinate (p) of a particle as a function of some other 
vector a (say the initial value of p) and of t. 

We first require the connection between a V and p V. We shall drop the affix a 
and retain p, so that now not Q(Si,Pi) but Q{V X , a{) 

= QU,0- 

Now S daV = S dp p V, 

but dp = -p\S daV\, 

S daV = -S da ViSpi p V, 

or since da is perfectly arbitrary 

V = -V Xl Sp lp V. (15) 

Thus operating upon equation (9) § 65 by V^piC ) and remembering that in that 
equation V must be changed to P V we get 

ViSpid- = V(v + P), (16) 

and this is our new equation of motion (Greenhill's (1) (2) (3) p. 448). 
In equation (1) § 12 above let us put for p + rj, p; and forp, a. 

Thus x^ = -ScoV .p. (17) 

Hence by equation (6d) § 14 above, 

strained vol. of el ./original vol. = S V 1 V 2 V 3 5'piP2P3/6, 

whence we see that 

DSV 1 V 2 V 3 Sp 1 p 2 p 3 = 6Do, (18) 

where Do is a constant which when a is taken as the initial value of p is the 
original density. This is the equation of continuity. 



"[Note added, 1892. At the time of writing the essay I did not notice that these equations are a 
particular case of the general equation for an elastic body already established (see equations (15m) 
§16and(15«)§16a) 

Dp =F - 2p[SVw<jw A.] 
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Cauchy's integrals of these equations. 

69. To obtain these integrals we require to express p V in terms of V. Now 
equation (15) of last section expresses V as a linear function of p V. The converse 
we have already seen how to get. In fact from equation (6h) § 3a above 

2/ p V = -Vp 1 p 2 SV 1 V 2 V, (19) 

where, with Greenhill, for brevity / is put for 5 , V 1 V 2 V 3 5'piP2P3/6. The only 
function we wish to apply this to is V p Vcr. We have 

2JV p Vo- = Vo- 3 VpipiS T1V2V3 

= (p 2 Spicr 3 -piSp 2 o- 3 )SV l V 2 V 3 , 

or interchanging the suffixes 1 and 2 in the last term 

JV p Va=p 2 Sp l a 3 SV 1 V 2 V 3 , (20) 

which gives the spin at any instant in terms of our present independent variables 
a and t. 

70. In order to obtain Cauchy's integral of equation (16) operate on it by 
VV( ). Thus 

= V(Vi + V 2 )Vi5 , pid- 2 = Spi& 2 . VV 2 V!. 

Now 

d(Spicr 2 . VV 2 Vi)/dt = So-icr 2 . VV 2 Vi + Spi& 2 . VV 2 Vi. 

The first term of this last expression is zero since the sign is changed by inter- 
changing the suffixes. From the last expression then 

Spicr 2 . W 2 Vi = constant 

= initial value of S£cr 2 W 2 £ = -W<x , 

where <r is the initial value of cr. Changing the suffixes and substituting in equa- 
tion (20) we have 

JV p Vo- = -p l SVo- V l , 

or giving / its value D /D from equation (18) and putting V p Vcr = e, W<x = e , 
we have 

e/D = -S(e /A))V.p. (21) 
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This is Greenhill's equations (4) (5) (6) p. 448. 

The physical interpretation of this equation is quite easy. Consider a small 
vector da drawn in the fluid initially. At the time t this will have become dp = 
-S daV . p. Thus we see that if a small vector ce /D be drawn in the fluid 
initially it will at the time t be ce/D, from which we infer that an element of a 
vortex filament will always remain an element of a vortex filament; or, a vortex 
filament or tube always remains a vortex filament or tube. Again we see that 
Te/D at any time varies directly with the elongation in the direction of the vortex 
filament so that Te varies as that elongation x the density, i.e. inversely as the 
cross-section of a small vortex tube at the point. This is not the easiest way of 
arriving at these results, but it is well to show in passing how easy of interpretation 
are our results. 

We see from equation (21) that if e = 0, e = 0. In other words, if the motion 
have a velocity potential at one instant it will have one always. 

Flow, circulation, vortex-motion. 

71. We are about to consider vortex-motion from another point of view, viz. 
that of circulation. 

In § 12 we saw that a strain due to a small displacement n could be decom- 
posed into a pure strain followed by a rotation, the vector rotation being Wn/2. 
If now for n we put the small vector <x 8t we see the propriety of calling V Vcr/2 
the spin. This therefore is taken as a definition of spin. Greenhill does not take 
this (usual) course but uses the property we shall immediately prove concerning 
circulation and spin to lead to his definition. 

It is not necessary here to define flow and circulation. Putting <x = V0 we 
see that for irrotational motion the flow = - J S crdp = f dcf) from one point to 
another is the increment in the velocity potential. Thus for mutually reconcilable 
paths it is always the same. 

Taking the circulation round a closed curve in the general case, the curve not 
inclosing any singular region of the fluid, we may transform the line integral into 
a surface integral by equation (8) § 6 above. Thus 

- / S cr dp = - ff S dX Vcr, (22) 

so that the circulation round the curve equals twice the surface integral of the 
spin. Hence Greenhill's definition of the spin. 
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72. From equation (9) § 65 we see that 

d ( r B \ r B r B 

I Scrdp] = - I Scr dcr - I S&dp 

dt\J A J J A J A 

pB pB 

= -i d(cr 2 )- d(v + P), 
Ja Ja 

I f S (T dp) = -[v + P + a 2 /2f A . (23) 



d ' "* 



dt 

Therefore for a closed curve 

d(JS(Tdp)/dt = 0, 

so that the circulation round the curve, and therefore the corresponding surface 
integral of the spin remains always constant. Taking the curve round a small 
vortex tube we once more arrive at the propositions enunciated in § 70 about 
vortices. 

Thus at any point of a vortex tube the strength which is defined as the product 
of the spin into the cross-section is constant throughout all time. Also it is the 
same for all points of the tube, for by equation (9) § 6 we have 

ffS dZVcr = fffSV 2 (Tds = 

for any portion of the tube. But the only part contributing to the surface integral 
is the ends of the part of the tube considered, so that the strength at these two 
ends is the same, and is therefore constant for the whole tube. 

73. These propositions are proved in the paper by Hicks already referred to. 
There is yet a third method of proof which, like Hicks's, is derived directly from 
the equation (9) of motion. 

We have from equation (1) § 62 

^-V - V— = -ScrV . V + VSo-V. = VjSo-iV. ; (24)* 

dt dt 

therefore dVVa/dt - VV& = V(Vi5 a x V 2 . cr 2 ). 



"[Note added, 1892. When first giving this in the Mess, of Math. 1884 and when again putting 
it in the present essay, I was unaware of Prof. Tait's paper in Proc. R. S. E. 1869-70, p. 143, 
where in the present notation he has for an incompressible fluid, (1) d/dt = d/dt - 5<rV [given 
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Now by equation (9) § 65, VV& = 0, so that 

dVVcr/dt = Wkt 2 SV 2 (t u (25) 

d(p + xVVcr)/dt = cr + xVVcr + xVV l cr 2 SV 2 cr l . 

Now ctiS V 2 cr 2 Vi = Va 2 V x S V 2 a x + WiV 2 S a x cr 2 + VV 2 cr 2 S Vi<n, but 
Wi V 2 S <T\(T 2 = 0, for an interchange of suffixes leads to a change of sign. There- 
fore 

Wi<t 2 S V 2 (Ti = VVcrS Vex - S VcrV . cr. 

Put now x = c/D where c is some small constant. Thus 

x/c = -DID 1 = -SVcr/D, 

and we get dip + cVVcr . /D)dt = cr - cS V<r V . cr . /D. (26) 

This shews that p' = p + cWcr/D is the variable vector of a particle, for the 
equation asserts that the time-flux of p' = the velocity at p'. Thus again we get 
the laws of vortices given in § 70. 

Irrotational Motion. 

74. We use this heading merely to connect what follows with what Greenhill 
has under the same on p. 449. It is not very appropriate here. 
For irrotational motion we have seen that we may put 

cr = V0. (27) 

If the fluid be incompressible we have further 

= SVcr = V 2 (f>. (28) 

Let T be the kinetic energy of this liquid. Thus 

-2T = DfffiVcf)) 2 ds = Dff<f>S dl V0 - Dfff(f)V 2 (f) ds 



in Cartesian form], (2) & = -Vv - Vp/D, (3) S Vcr = 0, (4) Vcr - d(Vcr)/dt = -SVcrV . cr, (5) 
Wo- = 0, (6) dVcr/dt = -SVcrV . cr. Perhaps I have not interpreted Prof. Tait's notation which is 
but briefly described correctly, but (4) should apparently be VV<x - d(Vcr)/dt = S VcrV . cr, and so 
agree with (5) and (6). It will be seen that the whole of §§ 512-13 of Tait's Quaternions, 3rd ed., 
is contained in this essay. I cannot at present recall whether this is owing wholly or in part to my 
being indebted to other old papers of Prof. Tait, or whether in writing his 3rd ed. he has arrived 
independently at the same treatment, but am inclined to the latter belief.] 
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by equation (9) § 6 above. Thus by equation (28) 

T = -±Dff<f>SdZV(f>. (29) 

Hence we see that if we have a vector cr which satisfies the equations V Vcr = 
and SVcr = throughout any singly-connected space, or more generally sat- 
isfies the second of these equations throughout any space, and also has zero cir- 
culation round any closed curve in that space*; and further satisfies the equation 
S dL cr = at the boundary; then the function 

fffcr 2 d* = 0, 

or cr = at all points of the space, for each element of the integral being essen- 
tially negative must be zero. 

75. The following important theorem now follows: If there are given; at 
every point of any region the convergence S V<x and the spin V V<x, at every point 
of the boundary the normal velocity (and therefore S <TLcr), and the circulation 
round every cavity which increases the cyclomatic number of the region; then the 
motion (as given by cr at every point) if possible, is unique. For let cr be one 
possible velocity, and if possible cr + r another. Thus t satisfies all the conditions 
that cr does at the end of last section, and therefore t = at every point or cr is 
unique. 

Motion of a solid through a liquid. 

76. We assume that there is no circulation of the liquid for any cycle; in 
other words, that if all the solids be brought to rest, so will be also the liquid. 

We shall take axes of reference fixed in the (one) moving solid. In the foot- 
note to § 67 is explained how the rotation of these axes is taken account of. The 
effect of the translation of the origin will be easily found. 

Let cr, co be the linear and angular velocities respectively of the moving solid. 
Let us put 

4> = -Serif/ - Scox, (30) 

where (f> is the required velocity potential of the liquid, and if/ and^f are two vector 
functions of the position of a point independent of cr and to. Let us see whether 



To insure that <p is single valued and so, that the part of the surface integral of equation (29) 
due to "barriers" is zero. 
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if/ and x can be found so as to satisfy these conditions, o~ and d are of course 
assumed as quite arbitrary. 

The conditions are first the equation (28) of continuity 

VV = 0, (28) 



which gives 



VV = 0,VV = 0, ( 3i) 



and second, the equality of the normal velocity -S U <iX V(p of the liquid at any 
point of the boundary with that of the boundary at the same point. Thus for a 
fixed boundary 

Scr{S dZ V . if/) + Sco(S dLV .%) = 0, 

whence on account of the arbitrariness of cr and to 

S dZV.if/ = S dZV.x = 0. (32) 

For a moving boundary again we have 

S dLV<p = S d~Z((T + Vtop), 
or Scr{dY, + S dY. V . iff) + Saj(VpdZ + S d~Z V .x) = 0, 

which gives 

dZ + S dLV.iff = 0, (33) 

VpdZ + S dZV.x = 0. (34) 

Now it is well known that if/ and x can be determined to satisfy all these 
conditions. In fact x being a coordinate of either iff or x these conditions amount 
to: — V 2 x = throughout the space, and S d~E Vx = given value at the boundary. 

77. We do not propose to find iff and^ in any particular case. We leave p. 454 
of Greenhill's article and go on to p. 455, i.e. we proceed to find the equations of 
motion of the solid. For this purpose we require the kinetic energy of the system. 
Calling this T we shall have 

T = -S(rZ(T/2-Su)<S>(T-Sa)Qco/2, (35) 
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where S, Q. and O are linear vector functions and the first two are self-conjugate. 
This is the most general* quadratic function of co and <x. It involves 21 indepen- 
dent constants, six in S, six in Q. and nine in O. When ift and x are known E, O 
and Q, are all known. We will obtain the expressions for them, and for simplicity 
will assume that the origin is at the centre of gravity of the moving solid. Let M 
be the mass of this last and p,co (where, as is well-known, ji is a self-conjugate 
linear vector function) the moment of momentum. Thus by equation (29) 

TT = -Dff(f>S dY. V(f> - Mcr 2 - Sojimd. 

Putting S cH,V(f) = S dL{cr + Vcop) at the moving boundary and zero at the fixed, 

TT = DfJ '(Scrip + SuxKSadi: + SojpdY) - Mo 2 -Scopoj 

where the surface integral must be taken only over the moving boundary. Thus 
noting that £ and D. are self-conjugate 

ZA = MA - \Dff(i//S Ad^ + d^S Aif/), (36) 

QA = fiA - \D§(xS Ap di: + Vpdi:S Ax), (37) 

-SA&A' = \Dff(S Ax S A' dZ + S ApdHS A'ifs). (38) 

These surface integrals can be simplified, for in each of these equations the 
first surface integral equals the second. In equation (36) in order to prove this 
we have merely to put for dL, -SdLV.il/ (equation (33)) when we shall find by 
equation (9) § 6 above (since we may now suppose the integrations to extend over 
the whole boundary), that 

ffif/SAdZ = -fffi// i SV l V 2 SAi// 2 d5 = ffdZSAif/. 
Similarly for equation (37). Again in equation (38), 

ffSAxSA' dH = -fffSV 1 V 2 SAxiSA'i{f 2 ds = ffSApdZSA'ifr. 
Thus finally for E, Q, and O 

ZA = MA- Dffi//S AdT. = MA- DffdZ S Aif/, (39) 

QA = fiA - DffxS Ap di: = ixA- DffVp JS S Ax, (40) 

-SA&A' = DffSAxSA'dH = DffSApdZSA'iJj. (41) 



*For let A(cr, <t)+B{<jj, <t)+C{oj, w)be this general function A, B, and C being scalar functions 
each linear in each of its constituents, and let 

So- = £{A(£, a-) + A(cr, £)}, (Dcr = £B(£, cr), Qw = ^{C(^, w) + C(w, £)}. 
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This last may be put in the following four forms 

®'A = -DffdZ SA X = -Dffif/S Ap dZ, (42) 

O/l = -DffxS AdZ = -DffVp dZ S At//. (43) 

Thus assuming that if/ and y are determined we have found T as a quadratic 
function of cr and to. 

78. If P, G be the linear and angular impulses of the system respectively our 
equations of motion are 

P + VcoP = F, 
G + VcoG = M, 

by the footnote to § 67 above. Here F and M are the external force and couple 
applied to the body. Now at the instant under consideration P = o-Vr, G = J7T. 
But if the origin had been at the point -p, P would still be ^VI whereas G would 
be JJT + VpP. Thus differentiation with regard to t does not affect the form of P 
but does that of G. In fact using the last stated values of P and G along with the 
last two equations and eventually putting p = cr, p = we get 

doVT/dt + Vco^VT = F, (44) 

dJJT/dt + Va)J7T + Vcr^VT = M. (45) 

Now from equation (35) we see that 

o-Vr = So- + O'w, (46) 

JVT = Oo- + Qa), (47) 

whence from equations (44) and (45) 

Ed- + O'cj + Voj Eo- + Vco O'w = F, (48) 

(Do- + Qco + Vcr Zcr + (Vcr <$>'u) + Vco dXr) + VcoQ.co = M. (49) 

Thus we see that with Quaternion notation even the general equations of mo- 
tion are not too complicated to write down conveniently. 

We now leave Greenhill's article and proceed to certain theorems not con- 
tained therein. The Cartesian treatment of these subjects will be found in Lamb's 
treatise on The Motion of Fluids, Chapters vi. and ix., which are headed Vortex 
Motion and Viscosity respectively. 
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The velocity in terms of the convergences and spins. 

79. We have seen in § 75 that if the spin and convergence are given for each 
point of a bounded fluid and the normal velocity at each point of the boundary, 
there is if any, but one possible motion. We shall see directly that this unique 
value always exists. 

At present we observe that we cannot find a motion giving an assigned spin 
and convergence for each point and any assigned velocity for each point of the 
boundary. If however with such data a motion be possible we can find the velocity 
at any point explicitly. By equation (19) § 10 above we have 

4ncr = -VfffVucrds, 

but by equation (9) § 6 

JJJVucr ds = JJu dLcr - JJJuVcr d$ 
4no- = ff[V]VudZo-- fff[V]VuVo-ds, (50) 

where the square brackets indicate that we may or may not retain the V at our 
convenience. This equation may be put 

4no- = ff[V](VudZo- - udZVcr) + fffuV 2 o-ds. (51) 

Both of these equations solve the question now proposed. 

If the fluid be considered infinite the surface integral of equation (50) vanishes 
if <x converges to zero at infinity and also that of equation (51) if in addition the 
spin and convergence converge to a quantity infinitely small compared with the 
reciprocal of the distance of the surface. 

80. We may now consider the case when Vcr (spin and convergence) is given 
at all points and S dLcr (normal velocity) at the boundary. It must be observed 
that the given value of the spin must be distributed in a solenoidal manner for 

SVVVcr = 0. 

Whenever the quaternion 

Anq = -JJjfVw Vcr ds = Vjjjfw Vcr ds, (52) 

is a vector, we see that all the conditions except the boundary ones are satisfied 
by putting a = q. Let us then see when q reduces to a vector. By equation (9) § 6 
above 

4nSq = -ffuS dUVcr + fffuS V 2 o~ds. 
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Thus q is a vector if this surface integral vanish. The surface integral vanishes 
if all the vortices form closed curves within the space. If they do not we must 
extend the space and make them form closed curves outside the original space. 
Extending the volume integral accordingly, we may put 

Ana' = -JffVuVads, (53) 

and now Vcr' = V<x. Suppose then that 

a = a' + a". (54) 

We thus get Vcr" = and may therefore put 

a" = V0, (56) 

where (p satisfies the equations 

V 2 = 0, (57) 

and S dLV(p = S dL(a - a') = known quantity. (58) 

Now it is well known that <p can be determined so as to satisfy these two 
equations. Therefore the problem under discussion always admits of solution. 

The above is equivalent to Lamb's §§ 128-131. His § 130 is the natural 
interpretation of the equation 

Ana = -fffVVuVads. 

His § 132 is seen at once from equation (50) above. For in the case he considers 
we, in accordance with § 7 above, take each side of the surface of discontinuity as 
a part of the boundary. Now [5 <TL a] a+h = so that for this part of the boundary 
we can leave out the part S <TLa and we get 

Ana = ffWu Vd^a- fffVVu Va ds, (59) 

so that if we regard [V dT,a] a +b as -2x an element of a vortex, we get the same 
law for these vortex sheets as for the vortices in the rest of the fluid. 



81. The velocity potential due to a single vortex filament of strength d6 is 
at once obtained by putting in equation (50) for Wads, 2d0dp. Thus calling a' 
the part of the velocity due to the filament 

2na' = -dOfVVudp = -deffu l VV l Vd'LV l 
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by equation (8) § 6 above. Now 

VV l VdYV l = -d^V\ + ViS rfSVx. 

Hence, since V 2 w = for all points not on the filament 

2na' = V(d8ffS rfS Vw). (60) 

Thus the velocity potential = (2ny l x the strength x the solid angle subtended by 
the filament at the point considered. 

Kinetic Energy. 

82. Let us put 

Anq = fffuVo-ds, (61) 

so that cr = Vq. (62) 

We may now find expressions for T, the kinetic energy, in one or two inter- 
esting forms. We have 

TT = -fffD(r 2 ds, 
whence IT = -fffDScrVqds = -ffDScrdXq + ^D x Scri V { qds. 

We assume that the fluid extends to infinity and that the vortices and conver- 
gences are all at a finite distance, so that at infinity cr is of order l/R 2 and q of 
order 1 /R. Thus the surface integral vanishes and we have 

2T = JffDiSo-iViqds, (63) 

or, putting in the value of q from equation (61), 

T = (87TT l fff 'fffuDiS cr l V l V 2 cr 2 ds x ds 2 . (64) 

These are Lamb's equations (28) and (29), p. 160, generalised. 
Similarly his equation 30 generalised is 

IT = JffDiSpa-i Vio-i ds, (65) 

for fffDiSpcr! Vio-i ds = fj DS per dY, cr - fffDS^cr^crds, 

by equation (9) § 6 above. But 

and the surface integral vanishes as before. 



§83.] HYDRODYNAMICS. 93 

Viscosity. 

83. To consider viscosity the assumption is made that the shearing stress 
which causes it is p. x the rate of shear of the moving fluid. // is assumed inde- 
pendent of the velocity and experiment seems to shew that it is also independent 
of the density. This last we assume though the variation with density can be easily 
treated. 

Consider a general strain x^- Since 

mSApv = Sx^XUXV = SA/Vx^XV*, 
Vxuxy = mx'^V/uv. 

Putting p and v for any two vectors perpendicular to to we see that the normal Uco 
to any interface becomes Ux'^to by the strain. Hence the interface co experiences 
a shear (strain) which equals the resolved part of xU to perpendicular to the vector 
X '~ l to equals resolved part of (x - x' l )Uco perpendicular to x'~ l(0 - Now when^ 
is the strain function due to a small displacement crdt, by § 12 above 

Xto = co - S co V . <x dt, (66) 

whence x'~ l( ° = to + V\Scoa\ dt, (67) 

and we see that the shear is the resolved part of 

-(S UcoV .cT + V^S Ucoctx) dt 

perpendicular to co, i.e. parallel to the interface to. From this we see by our as- 
sumption concerning viscosity that any element of the fluid is subject to a stress cp 
given by 

cpto = -Rco - n{ScoV . a + V\Scocr\) 

R being a scalar. Now we define the pressure p by putting 

3p = S£cf)£ = 3R + 2nS Vo- 
cpco = -pco + ^ficoS Vcr - fj(ScoV . tr + ViSco<Ti). (68) 

The equation of motion is 

D& = D¥ + tpA, 

or 

D(dtr/dt -ScrV .cr) = D& = D¥-Vp- fiVS Vcr/3 - //V V. (69) 

If n be not as stated independent of D this equation must be made to contain 
certain space fluxes of// which are quite easy to insert. 



"Kelland and Tait's Quaternions, chap. x. equation («). We have already used a particular 
case of this in § 15, above. 
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84. In § 179 of Lamb's treatise he considers the dissipation function due to 
the viscosity of a fluid in a way that seems to me misleading if not wrong. It 
appears as if he should add to the first expression of that section 

u dp xx /dx + v dp^/dx + w dp xz /dx. 

He refers to Stokes, Camb. Trans, vol. ix. p. 58. Let us give the quaternion treat- 
ment of Stokes's method. 

If T is the kinetic energy of any portion of the fluid 

Thus v d(Dd<5)/dt = 0, 

f = -JffDScr&ds, 
but D& = DF + <f>xV u 

so that T = -jJjDSo-Fds- jJJSo-frVidz 

or T = -fffDS(rFds-ffS(r0dS. + fffSo- l (f>V 1 ds. (70) 

If now <p be expressed in terms of p and p we have f depending on F, p and p. 
The part of t depending on F and p represents energy stored up as potential en- 
ergy of position and potential energy of strain respectively, but the part depending 
on p represents a loss of energy to the system we are considering the energy being 
converted into heat. 

Thus putting <f> = p + m we have by equation (68) 



(71) 



m(jL> = |yUd)5Vcr - p(StL>V . u + ^ \Six)(T\)\ 
= IfjuoS ffli; + 2fjuffa>, ) 

where if/ is given by equation (75) below. 

Thus we see that the rate of loss of energy is 

ffSvmdZ - fffS o-itjVi ds. (72) 

The surface integral is the work done by viscosity against the moving fluid at the 
boundary and the volume integral is considered due to the work done against the 
straining of the fluid. Thus we put 

F = -S(rnnV u (73) 

and call F the "dissipation function." 
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By equation (6) § 3 above 

F = -SMvt{, (74) 

where if/ is the rate of pure strain of the fluid, i.e. 

2ij/co = -ScoV .a -ViSaxTi. (75) 

Again by equation (18) § 18 above because F is quadratic in if/ 

F = -SM+aFa2, (76) 

so that from equation (74) we have, by § 4 above, 

m = ^QF/2. (77) 

Substituting for m from equation (71) in equation (74) 

F = -lfi(S(if f 2 -2M(K, (78) 

which gives F in terms of if/. 



Section VI. 
The Vortex- Atom Theory. 



85. If Quaternions can give valuable hints or indicate a promising method of 
dealing with the highly interesting mathematical theory of Vortex- Atoms, I think 
this alone ought to be sufficient defence of its claims to be within the range of 
practical methods of investigation. 

In what follows I think I may be said to have indicated a hopeful path to follow 
in order to test to some extent the soundness of this theory. 

Statement of Sir Wm. Thomson's and Prof. Hicks 's theories. 

86. Sir Wm. Thomson's theory is so well-known that it is not necessary to 
state it in detail. Matter is some differentiation of space which can vary its po- 
sition carrying with it so to speak certain phenomena, some of which admit of 
definite quantitative measurement. Perhaps the most important of these phenom- 
ena is what is called mass. Now, says in effect Sir Wm. Thomson, if we suppose 
all space filled with an incompressible perfect fluid the vortices in it are just such 
differentiations of space. They also carry about with them definite quantitative 
phenomena. Can we prove that these hypothetical vortices would act upon one 
another as do the atoms of matter, the laws of whose action are contained in the 
various Sciences, e.g. Physics, Chemistry and Physiology? The problem in its 
first stages at any rate is a mathematical one, but during the many years it has 
been before the mathematical world very little progress has been made with it. 

Hicks's extension of this theory is perhaps not so well known, but it seems 
to me quite as interesting and more likely to tally with the known phenomena of 
matter. He enunciates his theory in the Proceedings of the Cambridge Philosoph- 
ical Society, vol. in. p. 276. It differs from Thomson's simply in assuming that the 
fluid does not quite fill space — that there are in it bubbles*. These bubbles will 



*By "bubbles" of course I do not mean spaces occupied by another kind of fluid of smaller 
density but actual vacua. Thus a bubble may start into existence where none previously existed, 
or again a bubble may completely disappear. 
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find their way to where the pressure is least, i.e. speaking generally to the centre 
of some at least of the vortices. Thus we have another source of differentiation of 
space and general considerations seem to point to these differentiations being the 
true atoms, though of course "atom" is no longer a descriptive term. 

General considerations concerning these theories. 

87. In Maxwell's article Atom in the Encyc. Brit. p. 45, he says — "One of the 
first if not the very first desideratum in a complete theory of matter is to explain 

first mass and second gravitation In Thomson's theory the mass of bodies 

requires explanation. We have to explain the inertia of what is only a mode of 
motion and inertia is a property of matter, not of modes of motion. It is true that 
a vortex ring at any given instant has a definite momentum and a definite energy, 
but to shew that bodies built up of vortex rings would have such momentum and 
energy as we know them to have is in the present state of the theory a very difficult 
task. 

"It may seem hard to say of an infant theory that it is bound to explain gravi- 
tation." 

Now as Hicks tells us what induced him to give his theory was the promise it 
gave of explaining gravitation. But I believe nowhere does he point out the still 
more important result that probably on his theory we can explain inertia. 

The statement of the principal property of inertia put scientifically is that the 
motion of the centre of gravity of any two bodies approximates more and more 
nearly to uniform velocity in a straight line, the more nearly they are isolated 
from external influence. But this property is probably true of Hicks's bubbles. 
The centre of gravity of any portion of the fluid containing certain bubbles (1), 
(2) . . . (n) will if approximately isolated from all the rest of the fluid move approx- 
imately in a straight line, but this amounts to saying that the centre of volume of 
the n bubbles will also move uniformly in a straight line if the centre of the whole 
volume (bubbles and fluid) considered, move similarly. These conditions are not 
evidently true, but I think they are probably so when we consider groups of large 
numbers of bubbles. 

Whether this theory explains gravitation is one of the principal questions to 
be considered in the following first trial. 

Description of the method here adopted. 

88. Finding it practically impossible to consider the real problem of a num- 
ber of bubbles in the liquid, I consider the fluid continuous and not incompress- 
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ible. Let us assume that 

D = tanh(/Vc), (1) 

where c is some very small constant which in the limit = 0. For ordinary values 
of p, D very nearly = 1 . It is only when p becomes comparable with c that D 
varies. When p = 0, D = 0. Also 

P = J(dp/D) = c log sinhO/c)- (2) 

For ordinary values of p then, P = p but when p becomes comparable with c, 
P varies and when p = 0, P = -oo. If we assume then for the greater part of the 
fluid that p is large compared with the small quantity c we see that the fluid will 
have almost exactly the same properties as a liquid containing bubbles. We may 
now apply the equations of motion of § 64, § 65. 

Now we know the velocity at any point of an infinite fluid in terms of the spins 
and convergences at all points. From this we may deduce the acceleration in terms 
of the spins, convergences and time-fluxes. It so happens that by the equations 
in § 73 we can get rid of the time-fluxes of the spins. This greatly simplifies 
the further discussion of the problem. The equation that we thus deduce forms 
the starting point of our investigation, for the phenomena of gravitation, electro- 
magnetism, stress, &c, are exhibited and measured by the acceleration in bodies 
due to their relative positions. 

The equation we obtain at once gives a generalisation of the integral, equa- 
tion (12) § 65 above. 

At the end of this section I give an equation which is rather complicated but 
which promises to enable us to deal more rigorously with our problem than I 
profess to have done below. 

We proceed to the investigation just indicated. 

Acceleration in terms of the convergences, their time-fluxes, and the spins. 

89. Let us put for the convergence and twice the spin m and t respectively, 
so that 

SVo- = m, (3) 

Wo- = r. (4) 

We have seen that in equation (50) § 79 we may neglect the surface integral. 
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Thus 

47TO- = Vfffu(m + t) ds, (5) 

Anda/dt = Vfffud(m + r)/dt . ds. 
By §73 f = tSVct-StV .cr, 

but by equation (1) § 62 

dr/dt = t + S(TV . t, 

dT/dt = ScrV .t + tSV(T-StV . <r, 

or because S Vr = 

dr/d? = tS Act - aS At = V War (6) 

Andaldt = Vfffu(VVVo-T + dm/df) ds. (7) 

This equation is not in a convenient form for our purpose, for t and m may be 
and most probably will be discontinuous, so that it is advisable to get rid of their 
derivatives. Moreover it is advisable to allow the time-flux of m to appear only 
under the form d(m d%)/dt because in case of pulsations the time integral of this 
expression will be zero, whereas we can predicate no such thing of (dm/dt)ds, 
or indeed of m ds. Let us first then consider the first term in equation (7), viz. 
VfffuWVcrTds. Using equation (9) § 6, and neglecting the surface integral at 
infinity, this becomes 

- V JjjJ WuV(tt ds = VVVfffuVo-T ds [§ 9] 

= V 2 fffuV(TTds - VSVfffuVaTds 
= 4nVcrT + VfffScrTVu ds, 

by equation (19) § 10 above. 

Remembering now (§ 65 above), that 

& = dcr/dt - Vo-t - V(cr 2 )/2, 

we see by equation (7) that 

4n& = -2ttV . a 2 + VfffS ctt Vu ds + V fffu(dm/dt) ds, (8) 

whence by equation (9) § 65 above 

P + v-cr 2 /2 + (An^fffiScTTVu + u dm/dt) ds = H, (9) 

where H is a function of the time only. This is a generalisation of equation (12) 
§ 65, for assuming that t = we know by § 79 that dcf)/dt = (47r) _1 ffj u(dm / 'dt) ds. 
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90. This integral equation may be put into several different forms by means 
of equation (9) § 6 above. The form that is useful to us is the one in which instead 
of dm/dt we have d(mds)/dt as we have already seen. Now d(ds)/dt = -mds. 
Therefore 

d(m ds)/dt = (m - m 2 ) ds 

= (dm/dt - Sa Vm - m 2 ) rfs. 

Substituting from this for dm/dt, and then transforming by equation (9) § 6 
so as to get rid of the space variations of m involved in Scr Vm we get 

JJJu (dm/dt) ds = JJJu d(m d%)/dt + JJJ(um 2 - mSaVu - um 2 ) ds 
= JJJu d(m d%)/dt - JJJmS a Vw ds. 

Thus from equations (8) and (9) 

Ana = -2ttV . cr 2 + VfffSVu(Vcrr - ma) ds + Vfffu d(m ds)/dt, (10) 

and 

P + v - a 2 /2 + (An)' 1 fff{dsSVu(Var - ma) + ud(mds)/dt} = H. (11) 

Sir Wm. Thomson 's Theory. 

91. We are now in a position to examine the two theories. We first take 
Thomson's, which is considerably the simpler, and which therefore serves as an 
introduction to the other. We have then m = 0. Thus equations (10) and (11) 
become 

4na = -2nV . a 2 + V fffSarVuds, (12) 

p/D-a 2 /2 + (4ny l fffSaTVuds = H, (13) 

for in the present theories we may put v = 0. 

We shall consider the two terms on the right of equation (12) separately. The 
second term gives then an apparent force per unit mass due to a potential 

-(4ny l fffSaTVuds. (14) 

Comparing this with equation (33) § 46 above we see that this potential is the 
same as that of a magnetic system given by 

4ttI = VaT. (15) 
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Now this magnetism is zero where r is zero. In other words it is only present 
where the vortex-atoms are, and therefore it cannot be so distributed as to give an 
apparent force of gravitation, for taking the view of magnetic matter expressed 
in equation (34) § 46 we see that there must, in every complete vortex-atom, be 
a sum of magnetic matter exactly = 0. Nor again is it likely to explain the phe- 
nomena of permanent magnets, because, assuming that for any given small space 
including many vortex-atoms JjJVctt ds> is not zero, the apparent force produced 
will affect all other parts of space independently of whether this same integral 
for them is not or is zero. But to explain the phenomena of permanent magnets 
we must assume that the effect takes place only on portions of space where there 
is positive magnetic matter. This term then gives us no phenomena analogous 
to physical phenomena. As a matter of fact it probably has no visible effect on 
large groups of vortices, for there is no reason to suppose that the vector Vctt is 
distributed otherwise than at random. 

92. Let us now consider the other term in equation (12), and neglecting the 
term already considered, put 

& = -V.(T 2 /2. (16) 

The phenomena resulting from this are the same as would follow from a stress 
in a medium, the stress being an equal tension in all directions = -<x 2 /2. Now 
comparing equation (62) § 58 with equation (4) § 89, we see that cr depends on 
t/4tt in exactly the same way as does H on C. The question then arises — is 
the stress we are now considering equivalent in its mechanical effects upon the 
vortex-atoms to the stress given by equation (68) § 60 above, which explains the 
mechanical effect of one current on another? We saw in § 60 that this stress is 
a tension -H 2 /87r along the vector H and an equal pressure in all directions at 
right angles. The effects then would be the same only if <x be at right angles to 
the surface of our atom. But this is obviously not in general the case. From this 
analogy we can see however what approximately will happen to our atom. For 
instead of cr being at right angles to the surface it is in all probability very nearly 
tangential. Assuming that it is actually tangential we see that at the surface of 
the atom we have a tension exactly corresponding to the pressure which in the 
electric analogue will be exerted on this surface. In other words, the atoms will 
act on each other very approximately in what may be called a converse way to the 
small circuits in the electric analogue; i.e. where, in the electric analogue there is 
an attraction, in the hydrodynamic case there will be an apparent repulsion, and 
vice versa. 
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Now each vortex-atom forms a small circuit and therefore acts in the converse 
way to a small magnet. In other words, each atom acts upon each other atom as 
if it were charged with attracting magnetic matter. Thus we see that if we could 
suppose certain extra atomic vortices to exist and to be disposed throughout space 
in what at present must be considered quite an artificial manner with reference to 
the atomic vortices we could rear up a fabric which would explain gravitation. 
This conception however is of very little use for our present purpose. 

From these considerations I think we have every reason to believe that Thom- 
son's theory in its native simplicity does not promise to lead us to the physical 
phenomena of matter. We pass on therefore to Hicks 's. 

Prof. Hicks 's Theory. 

93. Hicks in his theory, as I understand him, assumes that the bubbles al- 
ways remain associated with the same particles of the fluid. This of course is 
probably not the case. By reason of the variation of the pressure with the time it 
is probable that evanescent bubbles start into existence and disappear at various 
parts of the fluid. This requires some few preliminary remarks. 

The particles of the fluid with which bubbles are permanently (i.e. throughout 
the greater part of each small but not infinitely small interval of time) associated 
are those where the intensity of spin is greatest. If the intensity of spin is quite 
various at different points we shall thus have vortices where there is generally no 
bubble, extra material vortices in fact. We must suppose these distributed quite 
at random till the more exact mathematical treatment of our problem leads us to 
suppose otherwise. Now evanescent bubbles will occur rather in these vortices 
than in parts of the fluid where is no vortex at all (if we may suppose such parts 
to exist), and of course they will occur more readily in stronger than in weaker 
vortices. At the present stage of the theory then we may suppose evanescent bub- 
bles to occur in all parts of the fluid. As a first approximation to the consideration 
of the effect of these bubbles we may assume a part m' of m to be continuously 
distributed through space. Putting 

m = m' + M, (17) 

we must suppose M to be present only at the material bubbles where it is probably 
discontinuous, whereas m' is continuous throughout both the material vortices 
and the rest of the fluid. 

Now when on account of variation of pressure m' and M are affected — is it 
probable that in the neighbourhood of a permanent bubble m' and M are of the 
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same or opposite signs? To answer this question observe that what we call m! 
continuously distributed is really a series of discontinuous values of m scattered 
at random through space so that m' is probably very small. A decrement of pres- 
sure will cause an increment of evanescent bubbles, i.e. a decrement in m' . An 
increment of the permanent bubbles will also take place the magnitude of which 
by what has been said concerning m' will not be by any means accounted for by 
the decrement in m' . There will therefore also be a decrement in M. Similarly for 
an increment in the pressure. M and m' may therefore be assumed to be of the 
same sign. 

After noticing that m' is continuous and therefore that there is no objection to 
introducing its space variations we are furnished with all the materials necessary 
for discussing our problem. The equation we shall use is (10) of § 90 above. We 
divide its discussion into two parts as follows. 

Consideration of all the terms except -V . (<x 2 )/2. 

94. The reasons that we have already seen in § 9 1 for neglecting Vfff S err Vw ds 
still hold good so we put this aside. This is not the case with -VjJJmScrVuds 
but we can neglect Vfffu d(m ds)/dt for the average value of d(m ds)/dt for any 
particle is zero. The only term to consider then is -VjJJmScrVuds. Putting as 
in equation (17) m = m' + M we have 

- V fffMS crVuds- Vfffm'S a Vw ds. 

The first term of this can be neglected for the same reasons as for neglecting 
that containing Vctt. Applying equation (9) § 6 to the last term and neglecting 
the surface integral as usual we get 

V fffuS V(mV) d* = Vfffu{m\m + M) + S(rVm}ds. (1 8) 

The last term can probably be neglected though we cannot give such good reasons 
as for the other terms we have neglected. At any rate in places not near permanent 
bubbles ScrVm' is as likely to be positive as negative and vice versa, so that 
such portions of space will on the whole produce no effect on the permanent 
bubbles. If Scr Vm' contributes anything for parts of space in the neighbourhood 
of permanent bubbles we must be content at present with the assumption either 
that the contribution is in general positive or that if it be negative it is not sufficient 
to cancel the effect of the positive term m'M. Remembering that m' is small 
compared with M we are left with the positive term m'M. This as can be easily 
seen from the form of equation (18) leads to an apparent law of gravitation for 
our permanent bubbles. 
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95. The gravitational mass which we must on this supposition assign to each 
permanent bubble varies as the average value of m'M for that bubble*. Now 
we saw in § 87 that the probable measure of mass of a permanent bubble was 
proportional to its average size. Do these two results agree? I cannot say, but 
even if they do not these considerations would still explain the motions of the 
solar system, but if the sun and Jupiter (say) were to collide their subsequent 
motion would not be that due to the collision of two bodies the ratio of whose 
masses is that which is accepted as the ratio of the sun's and Jupiter's. As a matter 
of fact however I should imagine that the average value of m'M for a permanent 
bubble is proportional to its average volume and this simply as a consequence of 
the reasoning in § 87 above. 

A conclusion at any rate to be drawn from the above is that there is a pre- 
sumption in favour of Hicks's theory explaining gravitation. 

Consideration of the term -V . (<x 2 )/2. 

96. In considering this term we adopt the method of § 92 and consider an 
electric analogue. The analogue is an electro-magnetic field for which in the 
notation of § 46 to § 60 above at every point, 

H = cr. (20) 

For this field we have at once 

4ttC = VVH = r. (21) 

The distribution of the magnetism in the field is somewhat arbitrary, but in the 
notation of equations (61) and (62) § 82 it will be found that everything is satisfied 
by putting 

Anl = -VSq. (22) 

This gives as it should S V(H + 4^1) = 0, which is in fact the only equation it is 
necessary to satisfy. We have further 

B = H + Anl = o- - VSq = VVq, 
so that A = Vq. (23) 



'There is one important difference to be noticed between this and Hicks's explanation of 
gravitation. His depends on the synchronous pulsations of distant vortices. I do not wish to 
imply that I do not believe in the existence of such synchronous pulsations, but by the above we 
see that gravitation can probably be explained independently of them. As a matter of fact such 
synchronous pulsations probably actually occur on account of the variation of H with the time. 
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Thus all the important vectors in the analogue are determined. It remains to 
compare the mechanical effects of the analogue with the term -V . <x 2 /2. 

I, it must be observed, is not confined to the bubbles, but is distributed 
throughout space. 

97. If we now assume that bubbles have not always existed in the positions 
which we call permanent, there cannot at the surface of the bubbles be any circu- 
lation round them. This makes the velocity at the surface almost normal to it, so 
that the stress given in equation (66) § 60 reduces to a tension -5'H(2B - H)/87r 
over the surface, i.e. we have a pressure 

-cr 2 /8n + S(rVVq/4n. 

Now on account of the absence of circulation Wq is very small and may therefore 
be neglected. Thus we get a pressure -cr 2 /8n, and are thus led once more to a 
"converse" of the analogue. This at once* leads to another reason for the law of 
gravitation if the pulsations are synchronous. This we have already seen to be 
probable. 

The present consideration of the subject is merely to point to a method of 
investigating the theory of vortex-atoms. I therefore leave the subject here, not 
attempting to force the phenomena we have been considering to tally with the 
known phenomena of electricity and magnetism. Nevertheless I may say that the 
prospect of discussing these things by means of the present subject can scarcely 
be considered as distant after what has gone before. 

To sum up, this first application of the method leads to a presumption in favour 
of Hicks's theory leading to an explanation of both the important properties of 
matter — inertia and the law of gravitation — and there is also reason from it to 
hope that the phenomena of electro-magnetism are not unlikely to receive an 
explanation. Thomson's theory on the other hand would seem to fail in the first 
two at any rate of those endeavours. 

98. We close the essay with the fulfilment of the promise made towards the 
end of § 88. In that section it will be remembered we considered a hypothetical 
fluid for which D = tanh(p/c), and made this do duty for a liquid containing 
bubbles. Strictly speaking our liquid is bounded at the bubbles and therefore as 
a bounded liquid should it be treated. For such a liquid we require an equation 



"[Note added, 1892. Because the density of attracting magnetic matter of the analogue = 
-SVH./4n = -ml An.] 
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corresponding to equation (10) § 90, and if possible equation (11) also. This last I 
have been unable to obtain, and I am not sure that to solve the problem explicitly 
is possible. 

Our problem only deals with an incompressible fluid, but as the removal of 
this restriction does not greatly complicate the work we will consider the general 
case of a bounded compressible fluid. We have 

Ana = V 2 fffuads = -VfffVuads 
= V fffu(r + m)ds- Vjfw dT. cr, 

whence 

Anda/dt = Vfffd{u(r + m)ds}/dt - Vffd(udLa)/dt. 

The justification of using d/dt on the left and d/dt under the integral sign will ap- 
pear if the increment (dcr/dt) dt in a at a given point in the time dt be considered. 
It will be observed that the meaning here to be attached to u will be -So- Vu, as 
in the differentiation d/dt with regard to the time the origin of u is assumed to be 
fixed. 

I shall now merely indicate the method of procedure. By the method exhibited 
in § 89 we can prove that 

d(iiTds)/dt = -VVw Vords, - uoSrAds 

and that 

d(um ds)/dt = -mS oVud% + u d(m ds)/dt. 

From this we can deduce that 

An dcr/dt = W V fffuVcrrd?, + Vfff{-mScr Vuds + ud(mds)/dt} 
- VffucrS TdZ- V JJk d(dl, cr)/dt + Vff <iX crS cr Va, 

and from this again we get 

o- = Vw + Vw' + Vv, (24) 

where w and vv' are scalars and v a vector given by 

Anw = -2ncr 2 + fff{ds S Vu(Vot - mo) + u dim ds)/dt}, (25) 
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so that w is in fact the H - v - P of equation (11) § 90 

4nw' = ff(S d^crSaVu-u dS dZ a/dt), (26) 

Anv = ff(-uoS rdZ + V dXcrScrVu - udV dZ cr/dt). (27) 

This last equation may be put in what for our purposes is the more convenient 
form 

Anv = -ffiuo-STdi: + d(uVdL(r)ldt}. (28) 

Again it may be put in a form free from d/dt; for V dL & = because the sur- 
face is a free surface and d{dY)fdt = ViS dL<T\*, as can easily be proved by 
considerations similar to those in § 83. Thus 

Anv = ff(-uo$TdL + VdLo$<rVu - uVV l( rS dL<r{). (29) 

In the problem we have to discuss m = 0, so that w gives only terms which 
we discussed in considering Sir Wm. Thomson's theory. Observing that if in w 
we change mds into -S dLcr we get for the part of w containing m, w'\ we see 
that W only gives terms that we have virtually discussed under Hicks's theory. 
We have however entirely neglected v. Are we justified in this? In the first place 
we have seen that if the bubbles have not always been associated with those parts 
of the fluid with which they now are there is round every bubble absolutely no 
circulation. This shows that for any one bubble Jjv dLa = 0\ and therefore we 
are justified in neglecting the last term in equation (28). We are probably also 



"[Note added, 1892. This should be d{dl)ldt = -mcTL + V^ dJ.a-\, and therefore equa- 
tion (29) should be 

Anv = JJ{VdI.crScr Vm - u(oSTdI. -mVd 1 Lo-+ VV icrS dLcr\)}. 

This does not affect our present problem because m = in our case.] 

1 One way out of many of proving this is as follows. Divide the bubble up into a number 
of infinitely near sections by planes perpendicular to the unit vector a. For any one section 
J S dpcr = 0. Consider d~L to be the element of the surface cut off by the following four planes, 
(1) the plane of section considered, (2) the consecutive plane of section, (3) the two planes per- 
pendicular to dp, and through the extremities of the element dp. Thus if x be the distance between 
the two sections, dp = x~ l VadI., whence 

ffSadI.o- = 0, 

for the surface of the bubble between the two sections. But adding, we may suppose this integral 
to extend over the whole bubble. Thus SaffV dZ cr = for the whole bubble; therefore a being 
a quite arbitrary unit vector we have for the whole bubble JJV dl.cr = 0. 
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justified in neglecting the first term, for probably t is very nearly tangential to the 
surface and therefore SrdL = 0. 

[Note added, 1892. The whole of this last section is in rather a nebulous 
stage, and since writing it I have not had sufficient leisure to return to the matter. I 
hesitated whether to include it in the present issue. But since, notwithstanding the 
absence of any reliable results, it serves very well to illustrate how investigations 
are conducted by Quaternions, I have thought it worth publishing. 

Should anybody feel inclined to attempt to apply the method or an analogous 
one it is well to note that in the Phil. Mag. June, 1892, p. 490, I have given the 
more general result sought in this last section. As the result is not there proved 
I give one proof which seems instructive. It exhibits the great variety of suitable 
quaternion methods of dealing with physical questions. It furnishes incidentally 
a fourth quaternion proof of the properties of vortices. It also illustrates how 
special quaternion methods developed for use in one branch of Physics at once 
prove themselves useful in other branches. 

Adopting the notation and terminology of Phil. Trans. 1892, p. 686, §§ 5-7, 
let cr and r be taken as an intensity and flux respectively, t still being = V Vcr 
and therefore r' = V Vcr'. Thus cr' and t' are the actual velocity and double spin 
respectively and the equation of motion is 

cr' = -V'(v + P). 

Putting cr' = x' 1(T an d operating on the equation by x' 

cr+x'd(x" l )/dt.cr = -V(v + P). 

Now 

X'd(x'- l )/dt . cr = -x'x'-'cr = - X 'cr' = V x Sa\a' = V . a' 1 12. 

& + V . cr' 2 /2 = -V(v + P), (A) 

which can easily be deduced from or utilised to prove Lord Kelvin's theorem 
concerning "flow," equation (23) § 72 above. 

As d/dt is commutative with V, t = or t is an absolute constant for each 
element of matter. This being interpreted at once gives the well-known properties 
of vortices in their usual form. 

If in the equation 4n(v + P) = S . V 2 JJJu(v + P) ds we carry one V across 
the integral sign, get rid of its differentiations which affect u by equation (9) § 6 
above, and then do the same with the other V we get 

4n(v + P) = JJ{(v + P)S dZ Vu - uS dZ V(v + P)} + fffuV 2 (v + P) ds. 
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At surfaces of discontinuity in cr, v and P will both be continuous, so that in- 
stead of JJ(v + P)S (TL Vw we may write JJ (v + P)S cTL Vw. In the last equation 
substitute throughout for V(v + P) from equation (A). Thus 

4tt(v + P) = ff h (v + P)S dZ Vu + ffuS d£ & - fffuS V&ds 

+ ( ffuS dH V . a' 2 12 - fffuV 2 cr' 2 ds/l) 
= ff h (v + P)S dZ Vu + ffuS dH&- fffuS Vcr ds + fffS Vu V . {a' 2 ) ds/2, 

which is equation (36) of the Phil. Mag. paper. 

If the standard position and present position of matter coincide it is quite easy 
to prove that 

dS dL' cr'/dt = S dl(& + yVjcrcri) 
d(SV(r'ds,')/dt = SV(&+VV l crcr 1 )d5. 

Substituting for S dl,&, S V<t from these in the last equation we get equation (32) 
of the Phil. Mag. paper; but this equation can also be proved directly. It should 
be noticed that equations (34) and (35) of that paper have been wrongly written 
down from equation (32). In each read + V . (<r 2 /2) for - V . (cr 2 /2).] 
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